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PREFACE 


This Final Report is submitted to the National Aeronautics and Space 
Administration's George C. Marshall Space Flight Center, Huntsville, 
Alabama, in response to the contract provisions of deliverable items 
associated with Structural Dynamics Payload Loads Estimates, Contract 
Number NAS8-33556. 

The study cook place during the period from August 1979 to October 
1982 under the direction of Mr. W. Holland of MSFC, Huntsville, Alabama. 

During this three year period the following documents were produced: 

1. Methodology Assessment Report 
August 1980, MCR-80-553 

2. Methodology Development Report 
August 1981, MCR-81-602 

3. Final Report 

September 1982, MCR-82- 601 

4. User Guide 

September 1982, MCR-82- 602 

5. Monthly Progress Reports 

The present Final Report together with the User Guide are intended to 
be largely self-contained. Chapter I deals with an overview of existing 
approaches to the problem of loads calculation. Chapter II introduces a 
full-scale version of an alternate numerical integration technique to 
solve the response part of a load cycle. A set of short-cut versions of 
this algorithm is developed in Chapter III. The implementation of these 
techniques is described in Chapter IV which deals with the software 
package. It should be noted that a complete documentation of the 
software can be found in the User Guide MCR-82-602 * Finally, in 
Chapter V the reader can find some tentative conclusions and a few 
suggestions of a more general nature. This report concludes with a 
rather extensive list of references. 

The author wishes to thank D. Devers, H. Harcrow and G. Morosow for 
their constructive comments and for reviewing parts of the manuscript. 
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INTRODUCTION 


The design and subsequent development of today's aerospace structures 
require extensive analytical and experimental studies. These 
investigations are necessary to ensure complete confidence in the ability 
of the total system to perform its required functions. Over the years a 
sufficient body of technical know-how has been accumulated to allow for 
an adequate system analysis prior to launch. However, with the 
appearance of larger and more sophisticated space structures, many of the 
analysis techniques are now being stretched to the limit of their 
capabilities. Indeed, the problem at hand sometimes has to be adapted in 
order to fit the analysis technique. This is particularly evident when 
dealing with large space structures and their corresponding finite 
element models. 

An important part of a design effort is the prediction of loads in 
the members of the structural system. A complete load cycle includes 
such items as modeling, modal reduction, modal analysis, load 
calculations ,response analysis, etc. As mentioned above a load cycle for 
a relatively small structure is fairly well standardized. The first step 
consists of constructing a discrete model for the structure. The analyst 
has to worry about the phenomena he wants to represent, how he will model 
certain structural elements, whft methods he should use, etc. If the 
structure is small, it is safe to say that an accurate structural model 
can be derived at reasonable cost. The next step in a load cycle is the 
calculation of the response of the structure to a given external force 
environment. Usually, a modal analysis technique is used, yielding a set 
of decoupled modal equations which are easily solved. Modes and 
frequencies essentially represent the physics of the structure and as 
such reveal much useful information besides the response. The last step 
in the load cycle is the actual calculation of the loads using the fact 
that the load vector is the product of the stiffness matrix ar the 
displacement vector. When the structure is small, all modes can be used 
to obtain the displacement vector and no acccuracy problems occur. The 
analyst finally sends the maximum loads to the stress engineer who 
calculates stresses and strains and finds out if certain structural 
members are correctly designed so they can survive the external force 
environment. When corrections in the design are necessary they will 
affect the overall response and therefore another load cycle is necessary. 

The techniques used on small structures are well established and 
tested on numerous real life structures. Large space structures pose a 
new set of problems primarily related to size and cost. Booster vehicles 
with payloads, such as the space transportation system aside from being 
large, have their own peculiar set of difficulties. The United States 
utilizes a rather small family of launch vehicles (boosters) to support a 
varied spectrum of satellite and spacecraft (payloads) programs. These 
launch vehicles have been carefully designed to accomodate a wide range 
of payload configurations. In general, the payload interfaces with the 
launch vehicle at a limited subset of candidate structural "hard" points 
at the payload launch vehicle separation plane. The Shuttle Orbiter is 
the latest example and is unique in the sense that it is reusable. The 
shuttle orbiter/payload(s) system therefore is not only a large space 
structure but part of it * namely the orbiter, is the same structure in 
every flight configuration. This makes the shuttle a unique challenge 
from the structural analyst's point of view. 



It is important that any candidate payload be designed to withstand 
the load environment transmitted to the payload from within the shielded 
payload compartment* Such environments commonly originate from a static 
(steady state) vehicle acceleration, a transient or dynamic event such as 
rocket motor ignition, or an acoustical environment* Very often, it is 
the transient dynamic response behavior of the payload that constitutes 
payload design load profiles; hence, it is important that proper 
attention be given to the payload transient response characteristics as 
influencing major design decision^, As an example let us consider the 
launching of the Shuttle Orbiter carrying a payload such as the Space 
Telescope. Obviously, when the rocket engines are ignited the Shuttle 
Orbiter will experience reaction forces. These forces will be 
transmitted to the Space Telescope in the cargo bay through the interface 
(i.e., through the connection points between the Orbiter and the Space 
Telescope). The space telescope then, will undergo elastic 
displacements. The question is, will the space telescope be able to 
withstand those vibrations without being critically damaged. The owner 
of this expensive telescope does not want to find out the answer to this 
question by trial and error. Therefore, an analysis is necessary. This 
analysis constitutes the load cycle as previously described. In the 
following we wish to point out some of the particular problems that arise 
during a load cycle conducted on a space structure such as the Shuttle 
Orbiter/Payload(s) /System. 

The first step in the load cycle, i.e., the modeling of the structure 
already poses a few unique problems. Even though research and 
development in the area of modeling techniques should continue, the 
engineer has available to him an excellent choice ot approaches. When 
the structure is large, the engineer has to take into consideration the 
size of the model. The question then becomes not so much do we have the 
necessary tools to model the structure? but: can we develop an adequate 
model of reasonable size and cost? To answer this question the engineer 
has to look at the particular phenomena he wants to describe and model 
accordingly. That is why the model developed by the structural engineer 
is often radically different from the model developed by the thermo 
dynamicist or the control engineer. Furthermore, if the booster model is 
reusable, the analyst must find ways to incorporate this established 
knowledge. Coordination between payload and booster organizations 
becomes critical. Reduction techniques must be considered. Is finite 
element modeling always recommended? Often, modeling cannot be done 
indiscriminately as not to result in oversized models, necessitating 
reductions of one kind or another. The analyst has to make sure he does 
not loose sight of the essential physical characteristics of the system 
in view of the large model. For example, is it necessary to model 
certain elements with two hundred degrees of freedom when the essential 
structural parameters only count five? Is there any benefit in hybrid 
models, partially discreet, partially continuous? 
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The next issue is the calculation of the response. Solving an eigen- 
value problem for the full-up discreet model becomes virtually impossible 
because of the enormous size of the model. The important question of 
model reduction becomes inevitable. Reducing a model can be done in 
several ways. A widely used technique is static reduction or Guyan 
reduction. This technique has to be administered with great care. Error 
analyses are scarce and much engineering judgement is necessary. Another 
popular technique is modal synthesis. Many different approaches and 
improvements of these approaches are available and will be discussed in 
more detail in subsequent sections. Again, no real error analysis is 
available and care must be taken not to loose vital information about the 
system. 

Once the modes and frequencies of the structure are known, the 
response can be calculated for several flight events. Many accurate and 
cost effective numerical integration techniques for the solution of a 
decoupled set of modal equations exist and have been used successfully. 
Large space structures should not present any special problem when the 
models can be reduced to acceptable sizes. The Space Shuttle however, 
poses a different kind of challenge. Indeed, the reusability of the 
shuttle makes it necessary to be able to deliver many load profiles at 
relatively short time intervals. The multitude of events, flights and 
payloads together with the reusability of certain models, flight data and 
other shuttle related knowledge, seems to require another breed of 
response techniques, indeed another approach to the entire load cycle. 
Short-cut methods in the load cycle which take complete advantage of the 
special nature of the shuttle problem should be developed. Much has been 
done but many possibilities still need to be looked at. Reduction 
techniques, response techniques, load calculation techniques, etc. should 
be investigated in the light of the particular needs of the shuttle. The 
relationship and interfacing between payload designer and payload 
integrator and booster organization and many other agencies and 
organizations is an integral part of this effort. Possibly, cookbook 
type approaches are necessary in order to facilitate the task of 
potential payload contractors or payload integration organizations. It 
seems appropriate to have short-cut techniques available in the 
preliminary design phase of a payload which are developed for specific 
circumstances and leave the payload designer with as much independence as 
possible. Only during the final verification cycle should expensive 
analysis tools be used. 

One of the major objectives of this report is to provide the reader 
with some insight into the multitude of problems and corresponding 
attempts at solving them which are available in the literature. The 
other objective is to develop a short-cut method in order to numerically 
integrate the equations of motion of the booster/payload system. 

Finally, we - 'sh to indicate also a few other possibilities for research 
in this quic j expanding field of payload integration techniques. 
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Finally, the actual calculation of the loads is also affected by the 
reduction of the model. In particular, when modal synthesis is employed, 
it is generally not possible anymore to use tie straight forward approach 
where the loads are given by the product of stiffness and displacement. 

It was found that significant error can result from such a displacement 
technique. Therefore, the so called acceleration approach is 
introduced. It makes the numerical problem more complicated but it 
yields accurate results. 


vii 



CHAPTER I. OVERVIEW OF EXISTING TECHNIQUES 
A. FULL-SCALE TEC.-HQUE S 


I* Introduction 

The objective of this chapter is to identify and discuss some of the 
most prominent techniques used in the course of a typical structural load 
cycle. This will acquaint the reader with state-of-the-art methodologies 
and the necessary background information in terms of a unified 
nomenclature. Also, many of the features of these methods can be in- 
corporated. into possible short-cut approaches. In addition, this 
chapter will give us the opportunity to more clearly identify the 
requirements of an acceptable short-cut methodology in connection with 
the STS (Space Transportation System). 

Any structural analysis starts with the derivation of the equatiomof 
motion for the system at hand. In the next section we shall derive a set 
of coupled differential equations describing the motion of a 
booster/payload system. 

2 . The Equations of Motion in the Discrete Time Domain 

The objective of this section is the derivation of the equations of 
motion - the booster/payload system. This set of equations subsequently 
must be solved in order to generate displacement histories which are 
needed for load calculations. 


Figure 1 shows the free body diagrams of the booster B and the 
payload P. The booster and the payload are connected to each other 
through the interface. Physically, the interface is the collection of 
structural "hard" points which the booster and th^ payload have in 
common. Mathematically, this means tMt the generalized interface 


displacement vector | j on the booster side of the interface must be 
equal to its equivalent j | on the payload side. Hence: 


K } * K I 


for all times t 


( 1 ) 
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Figure 1. Free Body Diagrams of 

Booster B and Payload P 
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Similarly, the generalized reaction vectors |r^ > and |r^ ? at the 
interface satisfy, ' ' 




for all times t 


( 2 ) 


From the free body diagrams in Figure 1 we can easily write the 
equations of motion for the booster B and the payload P as, 



(3) 


where {x«} represents the generalized displacement vector of the 
booster B. Furthermore, the matrices [Mg] and [Kg] are the mass and 
stiffness matrices of the booster B, respectively. Finally, the vector 
|Fg| represents the externally applied forces on the booster B and the 
vector |Rg|denotes the reaction of payload P on the booster B. Similar 
quantities for the payload are defined. 

Note that the top partition in equation (3) corresponding to the 
booster B, is still completely uncoupled from the bottom partition 
corresponding to the payload. In order to derive the equations of mot* 
for the coupled booster/payload system, we need to eliminate the a pi it 
unknown reaction vectors |Rg|and |Rp}« We shall now establish a 
convenient and physically meaningful way to accomplish this elimination. 
To this end, let us first introduce the interface/non-interface 
partitioning. Indeed, [Mg], [Kg], ^ Fq J and {Rb[ c n be 

partitioned according to interface and non-interface degrees of freedom, 
as follows: 
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Again, similar quantities can be written for P. Note that the subscript 
stands for Non-interface and I for Interface. Because there are no 
non-interface reactions we can write: 

K) - {») . K} ■ {°} <« 


substituting equations (4), (5) and (6) into equation (3) yields 



Both equations (3) and (7) represent the uncoupled equations cf motion of 
the undamped booster B and the undamped payload P. 

Next, let us solve the third partition of equation (7) for the 
non-interface displacement vector ( P \of the payload P. 

\*H f 

{<}■- [Cf [«a K}< * Kf ({<} - [C] K }- KJ H}) <»> 
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It Is now noted that the non-interface displacement vector{xjjjconaiste 
of two parts. To understand the physical meaning of these two terms let 
us assume that the interface displacements are aero i.e. £xf J - {oj . 

In that case it follows from Eq. '8) that the second term on the 
right-hand side can be Interpreted as the non-interface displacement of 
the payload with respect to the Interface. Let us denote this term by 
{xfj. It is then clear that the first term on the right-hand side of 
Eq. (8) represents the non-interface displacement of the payload due to 
the displacement {xj} of the interface only. Therefore, Eq. (8) can be 
written as : 


K} ■ M {<} ♦ ft) 

with 

[ s p] - - K» ] 1 Ki] 


It should be clearly uuderstood that < > is not the non-interlace 

displacement vector of the payload with respect to a fixed Interface Lut 
with respect to the interface (i.e., as seen by an observer moving with 
the interface). 

Equations (1), (9) and (10) can now be used to construct the 
following transformation: 



T3 NP 


oT j ! o 

”f"f" 

-°4M- 



^ 3 1 * + h 

of poo,; qomjJty 


(ii) 
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Note that we indicated the row and column sizes of the matrix partitions 
in the above equation. NB * number of booster non-interface dofs. NP * 
the number of payload non-interface dofs. IP * IB * number of interface 
dofs. 


Transformation (11) will eliminate the redundant set of interface 

process will also eliminate 
let us introduce a more 


Transformation (11) will eliminate the red 
displacements £ x j} in Equation (7) and in the 
the unknown reactions /R®l a ndiR P \ * First, 
convenient notation: I if Ilf 


NB 

no | o 1 


" 1 0 i 

ZB 

— > ” i 

Oil ! o 


0 xl 


— ! J 

i ; 


— +. 

NP 

0 1 S' I 


t 

T ! 

IP 

— J- — j 

01 I jo 


P ! 
1 


NB IB NP 




0 

0 

*p 

A '° 


& 

A 


( 12 ) 


with B = NB + IB, P * NP + IP and, 


OF POOR QOk 


M 

Px0 


°l s p 

— i — -I 

0| I 


PxnP 


f 

L oJzp 
NP 


(13) 


With this notation we now substitute the transformation (11) into 
equation (7) and premultiply by A T (T * transpose). This yields the 
following result: 


a 

p 


I x 

! W P 

NP 


Wi 



I I T 
K + T^iCT O I 
B PFP P P I 


T 

Wp 



r ^ 


X B 

- 4 

- — \ 


_P / 

> | 

*N 

V J 



(14) 
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where 


MB 

K} - i 

Bxl IS 


f b + Jr + F p 

F i + s p f n + F i 


OwiGiitAL r : ." r:: 

OF POOR QUALITY 


■ { g p} ■ Kl 


MPXl 


We also used equations (2) and (4). 

At this point a few remarks are in order* First we shall show tluat 
matrix product [Tp]^ • [Kp] • [Ip] is always zero* From equations 
(5) nnd (13) we have: 


r i T r i r i ! °1 T K ^-; <il 

M M H ■ _j 4 


‘Jkk 


TP P 

S P *nn + k in 


Substituting the expression (10) for [ Sp] into the lower half of 
the right hand side of equation (16) yields: 

[ s p] [4n] + [ k in] = ~ [ ^i] [C] Kn] + [ k in] = [°] 


because 


l>i»] ■ K I" ([«r] 


Secondly, the triple matrix product 

m r#. , r» P - . P 


0 * 0 


T 

si - 1 

p » 


VSi 

p i p 
k in! k ii_ 


is symmetric 


oj s 

- J _ 

I 

0 • I 


- .1 

1 p 

0 i K 
\ IN 


will be zero for a statically determinate inte^ r ~e. The interface is 
called statically determinate when the number of interface degrees of 
freedom is equal to the number of rigid bcdv degrees of freedom of the 
structure at hand. Otherwise, the inter tace is called statically 
indeterminate. To show that [K^] [Sp] + [ K j j ] in Equation (17) is 
zero for a statically determinate interface, let us first state that the 
numerical values of the elements cf this matrix are independent of the 
dynamical state of the structure. More specifically, we could constrain 
the booster and consider a case of static equilibrium for the system 
under the action of a force ^F^Jand all other applied forces equal to 
zero. As far as the stiffness matrix [Kp] is concerned, there will be 
no change. Therefore without loss of generality we consider the 
equilibrium equation of P under the action of J J • 
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or using the partitioned form of [Kpj , we can write 


W Hi + [■*] i x 

Ki kl + kj k p 


From equation (19) we can s^lve for J x N ( : 


j-4j ■ -KJ 1 ki] k | ■ W |4j (2D 

where we used equation (10). Substituting equation (21) into equation 
(20) yields: 

(k»] ki + ki])kj ■ HI ( 2 2) 

At this point we should note that when the interface is statically 
determinate no stresses can be set up in P by the interface 
displacements . Indeed, for a statically determinate interface the 

matrix [Sp] becomes a rigid body transformation, transforming the 
interface displacements into equivalent rigid body displacements of the 
non-interface degrees of freedom of P. Because, in addition we assu d 
that no other forces are acting on P, it is clear that{R^} is zero in 
Equation (25), from which it follows that 

[ K I»] M + Kl ■ 0 (23) 

This completes the proof. 

Finally, we note that in most applications of interest to this 
contract, the externally applied forces{F§\ , and { F ?1 are 

nonexistent. For example, STS payloads will be enclosed in the cargo bay 
and will not be exposed to external forces. Therefore we write 


Introducing Equation (24) into equations (15) and taking into account 
that rTpKpO- [0] we can write the final form of Equation (14): 

"B : T pVp| T pVp 1 W * TpKpTp j 0 ] k| _ U’ 

‘pVp i 'pVp JHT 0 ]$fy U‘t°; 


Wr l-i.l L° 


■ 0 - 
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represents tt^e generalized displacement vector of the free booster B. 

The vector contains all non-interface displacements of the booster 

and {xi} represents the interface degrees of freedom* Furthermore, 

the vector represents the non-interface displacements of the payload 

P with respect to the interface. The matrices [Mg] and [ Mp ] 
represent the mass matrices of the booster and the payload respectively 
and [Kg] and [Kp] represent the stiffness matrices. The matrix 
[Tp] is a constraint modal matrix, characteristic for the payload P. 

In case of a statically determinate interface [Tp] represents a rigid 
body transformation. Finally, [Fn} is the externally applied force 
vector acting on the booster B. The matrices [Tp]^ [Mp] [Tp] and 
[Tp] T [Kp] [Tp] contribute only to the interface degrees of 
freedom as can be seen from Equation (17) and from 



(28) 


with 



T 

Ti ^ matrix [T p ] [M^] [T p ] essentially represents the payload 

T 

mass reduced down to the interface. Similarly (Tp) [ Kp ] [Tp] 
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represents the payload stiffness reduced down to the interface. Note 
that when the interface is statically determinate no stiffness is 
transferred (Equation (23)). When the interface is statically 
indeterminate there are what is commonly called "constraint modes" , 
i.e., the interface displacements { not only induce rigid body 
displacements in the payload but also elastic deformations. This causes 
the triple product [Tpj* [Kp] [Tp] to be different from zero. 

It should al;'0 be noted that damping is not included in Equation 
(25). However, modal damping can always be taken into account later. 
Finally, the matrices [Mg], [Kg], [Mp] and [Kp] must be 
constructed by the appropriate organizations. The booster and payload 
organization will usually be different contractors. If in addition the 
integrator is another organization it is clear that there is an 
interfacing problem. A fair amount of coordination is necessary to make 
the transfer of information between these three organizations optimal. 
Unfortunately, this coordination is very difficult to establish resulting 
in considerable time delays. Moreover, these costs and delays repeat 
themselves for every load cycle (i.e., every time a change is made in the 
booster or payload). A typical example is the development of the Viking 
Orbiter System . Upward of nine organizations were responsible for 
hardware or integration functions which directly affected the evaluation 
of dynamic transient loads. The number of interfaces between those 
organizations resulted in difficulties in arranging for the necessary 
analyses at each organization, in obtaining the necessary data, in 
establishing priorities, in establishing output requirements and in 
correctly transferring data between organizations. The time duration for 
one load cycle ranged from three to twelve months which depended on the 
number of events, forcing functions per event, and complexities of the 
analysis. Of course, if the booster already has its final design, many 
of these problems can be avoided. Theoretically, only one transfer of 
booster data to the payload organization would be necessary. This 
concludes this section on the equations of motion of the booster/payload 
system in the discrete time domain. This is a necessary first step in 
any load cycle. 
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3. The Solution of the Equations of Motion 


a. The discrete time domain 


As stated in the Introduction, the objective of this study is to 
determine design loads for the payload structure* These design loads are 
then used to calculate stresses and strains that would exist in the 
structural elements that make up the payload* The stresses and strains 
allow the designer to determine the correct physical and geometric 
properties of the elements (mass, stiffness, lengths, cross sections, 
etc.) so that the structure does not fail when subjected to the external 


forces | l. An element loads equation is written as 

•MM 


(29) 


< p l 

in which 1 L € j is the load vector of an individual element e of the 
payload P, [k e ] is the stiffness matrix of the element, and [T e ] is 
the geometric compatibility transformation. The vector {xp} is the 
time dependent displacement vector of the payload satisfying Equation 
(25). Consequently, in order to determine { in Equation (29) we need 
to solve Equation (25). 


The most straightforward approach to determine [xp} is to solve 
Equation (25) as a set of simultaneous second order differential 
equations* There are several well established response routines that 
handle such problems (Runge-Kutta, Newraark-Chan-Beta Numerical 
Integration, etc*). This direct approach has the advantage of simplicity 
and accuracy. However, the main drawback is the high computational cost 
due to the large number of degrees of freedom used to describe today's 
aerospace models* In this connection, as a counter argument, one could 
ask the question if such large models are really necessary. In any case, 
unless a drastic reduction of the model is possible, this direct 
integration approach is not favored* In addition, this approach does not 
yield much qualitative information about the system* As will be 
explained below a modal approach reveals more qualitative information. 


b. Modal Analysis 

In this section we shall discuss a technique commonly known as modal 
analysis. This approach will lead us to an alternate solution method for 
Equation (25) and we shall show that it has some definite advantages over 
the direct solution of the set of differential equations (25) as 
discussed in Section 3a. 

We start the process with the homogeneous set of equations extracted 
from the top row of Equation (3). 


MN-WW-W 


(30; 
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Associated with this equation is an eigenvalue problem, 


(- A M + M)M ■ {°} 


( 31 ) 


where the vector^ j represents an eigenvector (mode shape) and an 
eigenvalue (natural frequency). The solution of this eigenvalue problem 
essentially produces a linear transformation matrix [4^1 (modal 
transformation matrix) in which each column represents a mode shape of 
the booster B. The main property of this modal transformation is that in 
the new normal coordinate system > the equations of motion (30) 

become uncoupled, i.e., if we apply the transformation, 


hi ■ NM 


(32) 


to Equation (30), and premultiply by [‘f’g]^ we obtain 

M T N Mhi + hr MM hi -|°t 

where the coefficient matrices of and j are now diagonal, 

M T MM ■ N- h] T MM- h] 


(33) 


(34) 


where [4>g] was normalized with respect to [Mg]. Equation (33) can 
then be written as 


+ tMh 



(35) 


The obvious advantage of applying the modal transformation Equation 
(32) is that Equation (35) now represents a set of decoupled independent 
second order differential equations that are easily solved. The price to 
pay however, is the cost of the solution of eigenvalue problem (31). 

There are many well established eigenvalue problem "solvers" available 
(Jacobi, Rayleigh-Ritz , etc.) 

The next step is to consider the homogenous equation, 

h] T H M hh hr hi M hi = 1"! 


ORiu.; > 
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Note that, 


ORIGINAL P/\Gi IS 

OF P''^ OVALITY 


M 1 M [ Ip ] -[<-]• N H 




(37) 


So that Equation (36) becomes 


(38) 


In the same way as we did for Equation (30) we can introduce a modal 
transformation, 



(39) 


with 


m T K.K] - i'i -kj 


(40) 


1^-2] . ^-P 2 1 

where we wrote L^pJ instead of! Ito simplify the notation. The modal 
matrix 1 <t> N J has as columns the "cantilevered" mode shapes of the payload P 


and 


t 


u)p j has the natural frequencies squared of the cantilevered payload 


(i.e., fixed interface) on its diagonal. Using Equations (39-40), we can 
write Equation (38) as 


<Ui 
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(41) 


Let us now apply the following transformation to Equation (25), 
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and preraultiply by 
(34) & (41) we obtain, 
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. Taking into account Equations 


I + ♦JXVb 
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( 43 ) 


At this point another eigenvalue problem could be solved for the 
coupled system represented by Equation (43). This would finally lead to 
a set of decoupled equations for the system. In fact, one could solve 
just one system eigenvalue problem for the homogeneous equation given by 
the homogeneous part of Equation (23). Numerically this approach does 
not yield much advantage over a direct solution. For one, the solution 
of a large eigenvalue problem is very exp'nsive and essentially increases 
with the cube of the dimension. The raoda, approach has the advantage of 
uncoupling the system equations, which allows us to investigate several 
load cases with little adued cost (if the structure and its constraints 
do not change). It also has the advantage of revealing physical 
characteristics. After all, modes and frequencies essentially embody all 
physical information. But their is potentially a much greater advantage 
present in the use of modes and frequencies. It can reduce the number of 
modal equations to be solved. This is why equation (43) was derived in a 
stepwise fashion. 

c . Modal Synthesis Techniques 

From what we have seen so far we can conclude that for small systems 
we have the necessary techniques to model the structure and subsequently 
find the response. The response can be found by direct integration of 
the discrete equations of motion or by first solving an eigenvalue 
problem. However, when the structure becomes large and the model has a 
very large number of degrees of freedom there is a problem. Direct 
numerical integration or solving a large system eigenvalue problem 
becomes prohibitive. Therefore, remedies have to be found. In section 
3b we already layed the groundwork for a reduction technique called 
"modal synthesis 11 . Indeed we introduced the free modes of the booster 
and the cantilevered modes of the payload(s). These mode groups are then 
"synthesized 1 * so as to yield the final coupled set of equations (43). As 
pointed out, this procedure by itself does not reduce the size of the 
problem and actually requires the costly solution of several eigenvalue 
problems. However, in most practical applications there is a possibility 
of defining a so-called "cut-off frequency". In these cases a Fourier 
series expansion of the force vector [ F^j shows that the energy content 
of the high irequency components is small compared to the energy 
contained ; n uhe low frequency components. Practically, this means that 
the respouoe of the structure due to the high frequency content of [ F^J“ 
can often be neglected. In this connection it should be noted that it is 
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relatively difficult to excite the higher modes of the structure to any 
large extent, especially when ^F*} only contains a few elements (i.e., 
only a few application points). The idea then is to first truncate the 
booster modes and payload modes in Equations (32) and (39) according to 
that "cut-off frequency". No thorough investigation has been published 
concerning the approximation involved in such a mode truncation. It is 
largely a "practical" matter supported by some theoretical 
considerations and the fact that it works. Truncation of the modes on 
the substructure level reduces the size of the eigenvalue problem related 
to Equation (43), Once this reduced eigenvalue problem is solved another 
truncation is possible, this time on the system level. To make this 
clear let us write equation (43) as follows: 


[M] \ q\ + f K ] £ q} * {q} 


(44) 


with 


[M] 


m 
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{ q} = ' 


( % 


0 


(45) 


where we should stress that [^g] and [^ N l are already truncated so 
that the size of [Ml and [K] is significantly reduced. Now, we solve the 
eigenvalue problem associated with the following homogeneous equation: 


[M] {q} + [V] {q} = {o} 

leading to 

lit - ml?} 

- m , m T unf] ■ r-a 2 j 


(46) 

(47) 

(48) 
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and 


(^} + - tVi T f q} <49) 


where we introduced the customary modal damping. We also truncated the 
modes Wi again according to the present cut-off frequency. The final 
set of uncoupled Equations (49) is now solved with a numerical integrator 
such as a Runge-Kutta method. 

Note that there are several synthesis techniques available, the 
purpose always being the reduction of the system size. When working with 
large structural systems and their corresponding large mathematical 
models, reduction methods become very important. Ways to reduce systems 
are: modal analysis, Guyan reduction, component mode synthesis, etc. 

The component mode synthesis approach as described in Section 3c is 
approximate in nature. Several techniques are available to improve on 
the accuracy, two of which we shall now discuss. 

d. The Residual Mass and Stiffness Method 


As explained in Section 3c in many cases it is possible to define a 
cut-off frequency which enables us to truncate the higher modes in 
Equations (32), (39) and (49) thereby reducing the size of this 
equation. Obviously, some accuracy in the response of the structure is 
lost due to the truncation of these higher modes. This loss of accuracy 
is especially apparent at the interface. The residual mass and stiffness 
method, instead of omitting these high modes will replace them with a set 
of "residual modes". The computation of these residual modes does not 
require any knowledge of the payload so that they represent a one-time 
computation effort not to be repeated as long as the booster stays the 
same. In order to determine the residual modes let us consider Equation 
(32): 


N 'N M 


(50) 


which represents the modal transformation for the booster B. Assuming a 
cut-off frequency was determined we can partition Equation (50) as 
follows : 



(51) 
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where l^g] represents the inodes with frequencies less than the cut-off 
frequency and those with higher frequencies. At this point one 

could neglect [<^g^]and calculate the response as a linear combination of 
the lower modes [<£ fi ]only. Usually this yields a poor accuracy in the 
response and the loads* The reason is that in most practical cases a 
significant part of the interface response is produced by the higher 
modes* Indeed, a typical interface is rather stiff and has little mass, 
i.e., that locally the interface has a high frequency content so that it 
responds significantly in the high frequency range. In truncating the 
higher modes the model does not include an adequate representation of 
that interface. The residual mass and stiffness method now, proposes to 
retain the static contribution to the response cf ■’brae higher modes. 

This leads to a much better representation of th'a inter ^ • The static 

contribution can be obtained from the following static C ion 

[■■] M 

derived from Equation (3). Substituting Equation (50) into Equation (52) 
premultiplying by [d^]T and recalling Equation (34) yields: 

[“» ] I’B j ’ [*b] 




Because we are only interested in the high frequency part, let us 
write Equation (43) as 



(54) 


So that from the bottom row in Equation (54) we have 


t-a m ■ w tij 


Finally, let us paitition [ i- ] in non-intertace and interface 
partitions, J 



(55) 


( 56 ) 
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Substituting Equation (56) into Equation (55) we obtain, 

[ 4 ] 


( 57 ) 


In principal we can use Equation (55) as it is and solve for 




B 

1 H 

p h 2 i r h i r 1 

P K 

R “ 

j q B| 

" L“b j k J - 
1 
1 


which can be substituted in Equation (51), yielding 




(58) 


(59) 


However, it should be noted that for every force component we keep, 
we add a degree-of-f reedom to the problem. If for example, £ contains 
many elements (i.e., many points of application) it may not pay off to 
use Equation (59), i.e., we may as well keep all the modes in T^uation 
(50). If however contains a small number of elements (for example, 

in case of a landing or a rocket motor ignition) we can use Equation (59) 
as it is, and obtain a better response for few added degrees of freedom. 
However, because the cut-off frequency was defined in such a way that all 
significant frequencies of are contained in the lower frequency 

range L, we can state that the booster model will adequately respond 
to and no significant portion of the response will be lost. 

Therefore, we can omit the term in £ in Equation (57) altogether and 
just keep the interface part in . The latter part in is 

important because will usually have a significant high frequency 

content (after all represents the effect of the payload on the 

booster and as such contains a wide range of frequencies). Because the 
interface usually has a high frequency content (as explained before) 
will induce a response at the interface primarily in the high frequency 
range which in turn will be, transmitted to the rest of the booster. 


On the other hand if contains reaction elements due to seme 

external constraints (e.g., a dock) we wish to retain these elements as 
well because they are eauivalent to elements of in the sense that 

they represent the unknown effects of the constraints and also, the 
interface between the constraint and the booster usually has a high 
frequency content (e.g., con actions between booster rna dock). 
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Ignoring the term in{F^| in Equation (57) and solving for yields 

K} ■ ] [+Bl] {'?} 


(60) 


which can be substituted into Equation (51), yielding 

{*4 •[♦»! [fltf’K] 1 ] 


-2 


(61) 


r, H ]M Mr t H 1 

The term ^ 9 hJ['h ^ I BlJ represent* the residual modes and they 

eplace [4>£] • Also, note that these modes only involve booster 
quantities which makes it a one-time computational effort. 

Let us now derive the modally coupled equations of motion for the 
booster/payload system* First, we substitute Equation (61) into the top 
row of Equation (3) and then we preraultiply by r * • - » -» « nT 

yielding 
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residual mass 


[♦a 


residual flexibility 


( 62 ) 


Before proceeding, let us consider the homogeneous equation extracted 
from the lower half of Equation (62) 

1 r .. tT ( ,.B 

i + l h ir,.H -| l« r n 

(63) 
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and solve the following eigenvalue problem 


leading to the modal t -ans format ion 
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K] ■ [♦.]{%' 


(64) 


(65) 
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and the properties 
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We shall make use of these properties in deriving the modally coupled 
equations of motion of the booster/payload system* To this end let us 
write Equation (25), 
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(67) 


Let us now introduce the fc Hewing notations 



so that, combining Equations (61), (65) and (68) we can write, 



(69) 


We now define a transformation similar to Equation (42) 
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where thi9 tirae(qg]and are given by Equation (68) 


Applying this transformation to Equation (67) and premultiplying 
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Using Equations ( 34) * 09) and (56) ve can show that 



so that Equation 171) becomes 



interface representation in the model wncre the model is subjected to a 
frequency cut-off. 


The result Equat ion (M) can then be used in any type of modal 
synthesis technique such as result Equation (74). Due to this interface 
improvement it is now possible to truncate the high booster modes while 
still obtaining an acceptable accuracy* In fact, inclusion of the 
residual mass and or stiffness into a modal synthesis technique seems to 
be the moat efficient approach currently available. Another verv good 
technique is discussed in the next section. 
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e. The Mass and Stiffness Loading Technique 

Another way of improving the interface representation in he booster 
model subject to frequency cut-off is given by a technique developed by 
Hruda and Benfieid and is based on Equation (25) which we repeat here for 
convenience , 
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Instead of solving eigenvalue problem (31) Hruda and Benfieid propose 
to solve the following eigenvalue problem. 


h * T fv P ] 4 |> 4 T Jvp])j*;J ■ j" 


(76) 


yielding the modal transformation 


w ■ W H 


\ 


(77) 


Again, the modal transformations (39) and (77) can be combined in 
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Substituting Equation (78) into Equation (75) and premult iply ing 
T y lelds , 
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where we used the properties 


‘i] T [ M B + T Jvp][*;] ■ [']• {*0 T [ k b + t p k p t p] [*i] 



Equation (79) now replaces Equation (43). The main difference lies 
in the fact that in solving eigenvalue problem (76), the booster 
interface is mass and stiffness loaded by [ TpMpTp ] and [TpKpTp] 
respectively; i.e., the booster interface is loaded with approximate 
dynamic effects^of the payload. In doing so, the new mode s [ 4 ^ ] and 
frequencies f (Oq Jwili inc 1 * le a good representation of the interface, and 
more resemble the system modes. This allow* us to reduce the number of 
booster modes in Equation (79) according tv the p^edeteraii * v,ut-off 
frequency. The disadvantage of this method in connection with the 
present study is that eigenvalue problem {lb) is dependent on the 
payload. This means that for every change in the payload we must solve 
this eigenvalue problem again although the booster does not change. This 
makes the Hruda/Benf ield technique less suitable for our purposes. 
However, if the changes in P are small, we can use the old booster modes 
as a first estimate to calculate the new booster modes in a Raleigh-Ritz 
type eigenvalue problem solver. 


f. The Coupled Base Motion Technique 


The coupled base motion technique as presented in this section is 
another way to obtain the response of a coupled booster/ payload system. 
This section will give us the opportunity to develop an alternative set 
of equations for equation (43). Indeed, we shall not only use 
"cantilevered" displacements for the payload P but also for the booster 
B, while only the interface will be free. The derivation is very similar 
to the one in Secton b. Let us define a transformation similar to 
Equation (9) but now for the booster B 
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Equation (11) can now be replaced by 
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Again, this transformation will eliminate the redundant set of 
displacements in Equation (3) and in the process it will also 
eliminate the unknown reactions and 

Introducing the notations 
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where this time, 
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Note that [I B j and [Ip] have different dimensions, namely (BXNB) and 
(PXNP) • 


Substituting transformation (83) into Equation (3) and premultiplying 
by A*** yields, 
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Equation (86) replaces Equation (25). 
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The basic idea for a base drive method is the separation of the 
booster response into two separate parts 
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part due to the action of <F >only 
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part due to the presence of the payload (•feedback) 
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ir is clear that vector (86) satisfies 
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The solution of Equation (90) is a one-time computational effort 
because it only involves booster quantities. If we now substitute 
Equation (87) into Equation (86) and take into account Equation (90) we 


obtain the following new set of equations 
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where and^x^ | are known from Equation (90) 

The second idea of a base drive method is to consider the bottom 
partition of Equation (91) and write it in the following form, 

[■JViJI’S, + ['JViJ i I --PpV P ](l ; fl + K"») 




If one is only interested in the design of the payload, Equations 
(92), (7) and (87) is all we need, to determine the response of the 
payload. If is known we can "base drive" the payload by the terms 

on the right-hand side of Equation (92) to obtain • Of 

course , is coupled into the booster equations in Equation (91). 

Equation (91) can also be written in terms of normal coordinates. To 
this end let us introduce the following transformation : 
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where anc * [^1 ] are °^ ta ^ ne< ^ f roni solving the following 

eigenvalue problems, 
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[^] T ^bVbH^] = tl3 . C4.® 1 T [Ib K B X b 3 C*S3 = t-S (96) 
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Substituting transformation (93) into Equation (91) and 
premult iplying by the transpose of the square transformation matrix in 
Equation (93) yields the modal form, 
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Equation (98) represents the coupled base drive equations of motion 
in modal coordinates* The main advantage of equations (98) lies in the 
force term of the righthanu side. If many different payload 
configurations must be investigated, then the force term is easier to 
evaluate than for example the force term in equation (43). This assumes 
of course that the booster and its interface remains the same so that the 
solution of Equation (90) remains the same. 


The main advantage of Equation (90) is that it is written in a form 
that can be used as a starting point for several short-cut methods. This 
will be discussed in the next chapter. 
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4. The Load Calculation 


As discussed in Section 3, the reason for solving the equations of 
motion of the booster/payload system is the determination of the 
displacement vector £xp 3 so that we can substitute this vector into 

I O - DsJCTjKl (99) 


in order to obtain the internal structural loads { } on an individual 
member e of the payload P. In principal Equation (99) could be used as 
is, but this "displacement " approach turns out to be very sensitive to 
inaccuracies in|xp^; e.g., truncating high frequency modes as we did 
in Section 3b could very easily lead to erratic loads 
Heurist ically speaking, {xp£ contains three parts, the static 
displacement, the rigid body displacement and the "vibrational" 
displacement • 


Therefore, if one has an error in { x P | one necessarily affect the 
accuracy of all three parts. For this reason one prefers the so called 
"acceleration method". Basically this approach is capable of separating 
the static and rigid body parts from the "vibrational" displacement. As 
a consequence one only makes errors in the "vibrational" part which often 
is the smallest part of the displacement vector £xp£- Such an 
"acceleration" approach which is consistent with modal synthesis 
techniques was developed by Hruda and Jones. [ 53 ] 


Recalling Equation (83) we can write 

{ X P } = [ T p] + 


(100) 



so that from Equation (99) we obtain 


<0 - OJ[tJ[t p ]{* i b | + [k ][I ][ij)2| 


From the bottom row of Equation ( 86 ) we obtain 
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and from the second row Equation ( 86 ) we obtain 
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Expressions (102) and (103) can now be substituted into Equation 
(101) yieloing an equation for in terms of accelerations. Many of 

the matrix multiplications involved in Equations (101-103) can be 
simplified by using a unit load solution, which is a feature of mr«>t 
finite element programs 

After this, it is possible to introduce cantilevered modes for the 
booster and the payload(s) and truncate these modes according to a preset 
cut-off frequency. In Chapter II we shall discuss the load calculation 
procedure in more detail. 
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5. Assessment 


There are very few publications dealing with the comparison of the 
different full-scale payload integration methods. In this section we 
shall discuss two such studies. A study by R. Hruda shovs that 
the residual mass and flexibility approach is one of the most effective 
in terms of cost, convenience and accuracy. The Craig/Bampton technique 
using interface payload(s) i;. also one of the best techniques. 

As a test structure, Hruda used two planar trusses coupled together 
by a statically indeterminate interface (Figure 2). Five different 
techniques were compared to the exact solution, i.e., the solution in the 
discrete time domain As discussed in Section 3a.: 

1. Hruda/Benf ield Technique - (Section 3e.): inertial coupling of 
truss-2 constrained modes onto free-free modes of truss-1 which 
was mass and stiffness loaded at its truss-2 interface 
degrees-of-freedom by the interface properties of truss-2. (IMSL) 

2. Craig/Bampton Technique (modal verson of Equation (86)): 
inertial coupling of truss-1 and truss-2 constrained modes onto 
a free-free modal representation of the interface 
degrees-of-freedom. (i/F) 

3. MacNeal Technique residual flexibility approach of coupling 
truss-2 constrained modes onto free-free modes of truss-1 which 
creates stiffness coupling (residual mass not included). (RFSWOM) 

4. Rubin Technique : (the residual mass and flexibility technique - 
secton 3d.): coupling of truss-2 constrained modes onto 
free-free modes of truss-1 which yields only inertial coupling, 
and, by consistent application to the mass and stiffness terms 
in the equations of motion, yields both residual stiffness and 
residual mass terms. (RFIWM) 

5. Rubin Technique but without residual mass const ribut ion for 
truss-1 . 

The truss problem as illustrated in Figure 2 represents a 
planar problem with three rigid body degrees of freedom (two 
translational, one rotational). Each pinned joint has two 
translational degrees of freedom. The interface is statically 
indeterminate becau' - there are six interface degrees of 
freedom. The heavy masses (asymmetric with respect to 
interface) are added to produce interface distortion. The 
forcing function is a ramp function. (RFIWOM) 

The "exact 11 results, against which all comparisons were made, 
were obtained by extracting eigenvalues , eigenvectors, and loads 
directly from a finite element discrete/physical model using no 
modes at all. 
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Five different cases we investigated 
EXACT: Discrete modal 70 DOF 

CASE A: ModcJlly coupled, 70 modes retained (*1002) 

CASE B: Modally coupled, 50 r- J es retained (*712) 

CASE C: Modally coupled, 32 modes retained (*462) 

CASE F: Modally coupled, 19 modes retained (=272) 

Hruda used the following comparison values: 

Frequencies: Percent error against the "exact" solution 

Modes: An error vector is formed^^- ) and its 

norm is calculated (which is defined as the Root 
Square Sum of the the elements of the vector); 
the comparison value is then defined as the norm 
of the error divided by the norm of the 
base/exact mode. Note that the norms are based 
on the modal amplitudes of all degrees of freedom 
from the coupled system. 

Loads: Loads were calculated at the truss interface on 

both the truss-1 and truss-2 joints. A percent 
error of the absolute value of the largest 
{’•'ther maximum or minimum) value from a given 
case against the absolute value of the largest 
value from the exact solution. 



where E * Exact, and N * Case being compared. 
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The results are presented in Tables 1-12 

For the 100% case-A, the MacNeal technique requires the inversion of the 
residual flexibility matrix to obtain a "residual stiffness"* When 
attempting to retain all (100%) of the modes, this residual flexibility 
matrix is a function of the interface highest frequency modal amplitudes 
which can cause an ill-conditicned matrix (as in the present case). 

Since this is an unrepresentative case, it should not be deduced that 
this is an unacceptable technique. As can be seen in succeeding cases, 
where more residual modes are available, the MacNeal technique falls into 
line with other techniques. Note that in cases B, C, and F, in both the 
frequency and mode shape comparisons, that the Macueai and the Rubin 
technique without residual mass are identical, thereby numerically 
supporting R. Coppolino's contention that these two techniques are 
equivalent for modal synthesis. In comparing the loads it is seen that 
the MacNeal column for case-A reflects the propagation of the 
ill-conditioning mentioned earlier. Loads were calculated by the modal 
acceleration technique (Section 4) for all methods except for the MacNeal 
technique. 

Due to the stiffness coupling involved in the MacNeal method, a 
complete modal acceleration technique for calculating loads could not be 
used, therefore, the modal displacement techniques of calculating loads 
was used. Because of this, the larger loads inaccuracies for this method 
must be attributed to the method of loads calculation and not to th'„ 
method itself. 

In conclusion ve can state that methods 1 through 5 are acceptable. 
However, the Rubin Technique (Residual Mass and Stiffner- Approach) seems 
to outweigh the other approaches in terms of cost ' A *nience. 

Again, it should be noted that this method does roc any knowledge 

of payload properties which makes it very valuable fer &nal>jis ol 
STS- appl icat ions . 



Figure 2 Structure Used for Comparing Coupling Techniques 


ORIGINAL PAGE H 
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X 

Table 1 
Case A 
Diff for 

: Frequency Cooper i son 

■ 100 Z of Available Modee 
Various Modal Coupling Tech 

OR.'GuN'Ac. i 
OF POOR QU, 

ali j y 

Mode 

No 

Exact Freq 
( Hr ) 

IMSL 

1 /F 

RFSU 0 M 

RFIWM RFIWOM 

4 

.76 

-.00 

-.00 

- 49.17 

-.00 

-.00 

5 

1.76 

.00 

-.00 

6.76 

-.00 

.00 

6 

2.84 

-.00 

-.00 

- 21.81 

-.00 

-.00 

7 

3.08 

-.00 

-.00 

- 3.63 

-.00 

,00 

8 

3.80 

.00 

-.00 

- 9.24 

- . 00 

.00 

9 

4.62 

- .00 

-.00 

-16.35 

-.00 

.00 

10 

5.11 

-.00 

-.00 

- 24.08 

-.00 

.00 

11 

5.50 

-.00 

-.00 

- 6.02 

-.00 

.00 

12 

5.81 

-.00 

-.00 

- 3.31 

-.00 

.00 

13 

7.69 

.00 

-.00 

- 6.08 

-.00 

.00 

14 

8.69 

-.00 

-.00 

-.60 

-.00 

.00 

15 

9.14 

-.00 

-.00 

- 3.32 

-.00 

.00 

16 

9.42 

-.00 

-.00 

- 1.76 

-.00 

.00 

17 

9.73 

-.00 

-.00 

-.92 

-.00 

,o0 

18 

9.85 

-.00 

-.00 

\ 

CO 

-.00 

-.00 

19 

10.36 

-.00 

-.00 

- 1.36 

-.00 

.00 

20 

10.43 

-.00 

-.00 

-.33 

-.00 

.00 

21 

10.79 

-.00 

-.00 

- 2.°9 

-.00 

.00 

22 

10.90 

-.00 

-.00 

- 2.55 

-.00 

.00 

23 

11.37 

-.00 

-.00 

- 1 .85 

-.00 

.00 

24 

11.49 

-.00 

-.00 

.80 

-.00 

.00 

25 

11.78 

-.00 

-.00 

.39 

- . 00 

.00 

26 

11.96 

-.00 

-.00 

-.12 

-.00 

.00 

27 

12.03 

-.00 

-.00 

.19 

-.00 

.00 

28 

12.20 

-.00 

-.00 

.40 

-.00 

.00 

29 

12.43 

-.00 

-.00 

.12 

-.00 

.00 

0 

12.50 

-.00 

-.00 

5.58 

-.00 

.00 

31 

12.75 

-.00 

-.00 

4.00 

-.00 

.00 

32 

13.29 

-.00 

-.00 

2.54 

-.00 

.00 

33 

13.51 

-.00 

-.00 

3.12 

-.00 

.00 

34 

14.22 

-.00 

-.00 

-.19 

-.00 

.00 

35 

14.53 

-.00 

-.00 

.31 

-.00 

.00 

36 

14.86 

-.00 

-.00 

-.11 

-.00 

.00 

37 

15.19 

-.00 

-.00 

1.59 

-.00 

.00 

38 

15 . ^4 

-.00 

-.00 

1.23 

-.00 

.00 

39 

15.69 

-.00 

-.00 

2 . 5 C 

-.00 

.00 

40 

16.16 

-.00 

-.00 

.55 

-.00 

.00 

41 

16.17 

-.00 

-.00 

.66 

-.00 

00 

42 

16.28 

-.00 

-.00 

3.34 

-.00 

.00 

43 

16.86 

-.00 

-.00 

.46 

-.00 

.00 

44 

1/.07 

-.00 

-.00 

.21 

-.00 

.00 

45 

17.17 

-.00 

-.00 

5.12 

.00 

.00 

46 

18.17 

- . 00 

-.00 

.03 

-.00 

.00 

47 

18.17 

-.oo 

-.00 

10 . «7 

.00 

.00 
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Table 1: Frequency Comparison (Cone 1) 

Case A * 100Z of Available Modes 
Z Diff for Various Modal Coupling Tech 


Mode 

No 

Exact Freq 
(Hz) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

48 

20.11 

-.00 

-.00 

.81 

-.00 

.00 

49 

20.27 

-.00 

-.00 

3.02 

-.00 

.00 

50 

21.11 

-.00 

-.00 

.14 

-.00 

.00 

51 

21.15 

-.00 

-.00 

.08 

-.00 

.00 

52 

21.26 

-.00 

-.00 

.07 

-.00 

.00 

53 

21.29 

-.00 

-.00 

.53 

-.00 

.00 

54 

21.45 

-.00 

-.00 

.02 

-.00 

.00 

55 

21.47 

-.00 

-.00 

2.91 

-.00 

.00 

56 

22.06 

-.00 

-.00 

.45 

-.00 

.00 

57 

22.15 

-.00 

-.00 

2.18 

-.00 

.00 

58 

22.63 

-.00 

-.00 

5.09 

-.00 

.00 
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Z 

Table 2 
Case B 
Diff for 

: Fr 

- 71 * 
Vario 

uency Comparison 
Available Modes 
Modal Coupling Tech 



Mode 

No 

Exact Freq 
(Ht) 

IMSL 


I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

.76 

.00 


.00 

.00 

.00 

.00 

5 

1.75 

.00 


.oc 

.00 

.00 

.00 

6 

2.84 

.00 


.00 

.00 

.00 

.02 

7 

3.08 

.00 


.00 

.00 

.00 

.06 

8 

3.80 

.00 


.00 

.00 

.00 

.00 

9 

4.62 

.00 


.00 

.04 

.00 

.04 

10 

5.11 

.00 


.00 

.00 

.00 

.00 

11 

5.50 

.00 


.00 

.05 

.00 

.05 

12 

5.81 

.04 


.05 

.03 

.03 

.03 

13 

7.60 

.00 


.00 

.21 

.00 

.21 

14 

8.69 

.02 


.02 

.03 

.01 

.03 

15 

9.14 

.04 


.06 

.12 

.04 

.12 

16 

9.42 

.01 


.01 

.01 

.00 

.01 

17 

9.73 

.01 


.02 

.02 

.01 

.02 

18 

9.85 

.03 


.04 

.03 

.03 

.03 

19 

10.36 

.01 


.01 

.06 

.01 

.06 

20 

10.43 

.00 


.00 

.01 

.00 

.01 

21 

10.79 

.03 


.05 

.04 

.02 

.04 

33 

10.90 

.03 


.03 

.04 

.02 

.04 

23 

11.37 

.00 


.00 

.03 

.00 

.03 

24 

11.49 

.00 


.01 

.05 

.00 

.05 

25 

11.78 

.03 


.03 

.08 

.03 

.08 

26 

11.96 

.08 


.09 

.08 

.07 

.08 

27 

12.03 

.02 


.03 

.04 

.02 

,04 

28 

12.20 

.13 


.13 

.10 

.10 

.10 

29 

12.43 

.01 


.01 

.02 

.no 

.02 

30 

12.50 

.03 


.09 

.07 

.03 

.07 

31 

12.75 

.07 


.07 

.08 

.06 

.08 

32 

13.29 

.04 


.05 

.06 

.02 

.06 

33 

13.51 

.oe 


.07 

.07 

.05 

.07 

34 

14.22 

.13 


.13 

.78 

.09 

.78 

35 

14.53 

.09 


.13 

.29 

.07 

.29 

36 

14.86 

.01 


.01 

1 17 

.01 

.17 

37 

15.19 

.06 


.07 

.11 

.05 

.11 

18 

15.54 

.01 


.01 

1.17 

.04 

1.17 

39 

15.69 

.02 


.02 

1.61 

.03 

1.61 

40 

1 C . 16 

.01 


.02 

.13 

.09 

.13 

41 

16.17 

.07 


.08 

2.68 

.12 

2.68 

42 

16.28 

.37 


.44 

3.88 

.34 

3.88 

43 

10 .86 

. 13 . 


.12 

7.73 

.25 

7.73 

44 

17.07 

.06 


. 10 

18.31 

i . 3 2 

18.51 

4 5 

17.17 

t .00 


1 .00 

66.94 

3 . 86 

66.94 

46 

18.17 

.02 


.02 

72.29 

9.54 

72.29 
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Table 2: Frequency Comparison (Cone l ) 




X 

Case B 1 
Diff for 

■ 71X of Available Modes 
Various Modal Coupling Tech 



Mode 

No 

Exact 

(Ht) 

Fraq 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4> 

18.17 


.03 

.04 

122.22 

11.51 

122.22 

48 

20.11 


.07 

.09 

473.05 

7.55 

473.05 

49 

20.27 


.31 

.3 

942.70 

8.65 

902,70 

50 

21.11 


4.07 

4.20 

536.67 

19.40 

536.07 
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Table 3: Frequency Comparison 

Case C ■ 46% of Available Modes 
% Diff for Various Modal Coupling Tech 


Mode 

No 

Exact Freq 
(Hx) 

IMSL 

l/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

• 76 

.00 

.00 

.02 

.00 

.02 

5 

1 *75 

.02 

.00 

.00 

.00 

.00 

6 

2.84 

.04 

.02 

.12 

.01 

.12 

7 

3.08 

.01 

.00 

.28 

.00 

.28 

8 

3.80 

.02 

.00 

.00 

.00 

.00 

9 

4.62 

.12 

.05 

.26 

.02 

.28 

10 

5.11 

.02 

.01 

.01 

,00 

.01 

11 

5.50 

.26 

.10 

.35 

.05 

.35 

12 

5.81 

.16 

.24 

.12 

.08 

,12 

13 

7.69 

.03 

.04 

1,55 

.06 

1.55 

14 

8.69 

.14 

.22 

.17 

.08 

.17 

15 

9.14 

.32 

.33 

.41 

.17 

.41 

16 

9.42 

,05 

.07 

.07 

.01 

.07 

17 

9.73 

.13 

.18 

.34 

.08 

. 34 

18 

9.85 

.30 

.27 

.31 

.18 

.31 

19 

10.36 

.11 

.11 

.51 

.12 

,51 

20 

10.43 

.09 

.04 

.45 

.02 

.45 

21 

10.79 

.17 

.33 

,73 

.11 

.73 

22 

10.90 

.41 

.42 

3.43 

.55 

3.43 

23 

11.37 

.06 

.04 

3.23 

1.49 

3.23 

24 

11.49 

.04 

.11 

5.14 

3.66 

5,14 

25 

11.78 

.17 

.17 

5.52 

2.76 

5,52 

26 

11.96 

.40 

.90 

4.79 

3.70 

4.79 

27 

12.03 

.30 

3.16 

53,65 

4.71 

53.65 

28 

12.20 

1.73 

2.35 

77.63 

3.91 

77.63 

*>9 

12.43 

.41 

4.66 

80.30 

14.58 

80.30 

30 

12.50 

1.94 

12.32 

143,04 

22.38 

143.04 

31 

12.75 

3.52 

14.90 

167.76 

29.91 

167.76 

32 

13.29 

11.88 

61.19 

630.91 

63.27 

630.19 
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Table 4: Frequency Comparison 

Case F * 27Z of Available Modes 




Z 

Diff for 

Various Modal 

Coupling Tech 



Mode 

No 

Exact Freq 
(Hz) 

IMSL 

I/F 

PFSWOM 

RFIWM 

RFIWOM 

4 

.76 


.00 

.00 

.05 

.00 

.05 

5 

1.75 


.02 

.00 

.02 

.00 

.02 

6 

2.84 


.06 

.06 

.31 

.03 

.31 

7 

3.08 


.07 

.00 

4.32 

.19 

4.32 

8 

3.80 


.03 

.01 

.03 

.00 

.03 

9 

4.62 


.18 

.14 

.71 

.08 

.71 

10 

5.11 


.05 

.03 

.03 

.01 

.03 

11 

5.50 


.46 

.34 

1.64 

.49 

1.64 

12 

5.81 


1.03 

1.13 

23.29 

4.70 

23.29 

13 

7.69 


.14 

.16 

15.39 

8.20 

15.39 

14 

8.69 


1.05 

1.14 

45.28 

1.30 

45.28 

15 

9.14 


1.20 

6.85 

84.48 

8.49 

84.48 

16 

9.42 


.37 

27.17 

98.00 

27.57 

98.00 

17 

9.73 


2.23 

36.67 

97.08 

36.41 

97.08 

18 

9.85 


18.25 

37.17 

203.74 

37.43 

203.74 

19 

10.36 


22.07 

103.62 

279.52 

105.23 

279.52 
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Table Si Mode Shape Comparison 
Case A • 100Z of Available Modes 



% 

Diff for 

Various Modal 

Coupling Tech, 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

l/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.00 

.00 

104.09 

.00 

.00 

5 

13.95 

.00 

.00 

92.76 

.00 

.01 

6 

13.45 

.00 

.00 

104.91 

.00 

.01 

7 

13.70 

.00 

.00 

94.49 

.00 

.02 

8 

12.23 

.00 

.00 

132.95 

.00 

.02 

9 

11.21 

.00 

.00 

100.73 

.00 

.03 

10 

9.54 

.00 

.00 

153.35 

.00 

.05 

11 

8.75 

.00 

.00 

104.39 

,00 

.12 

12 

13.25 

.00 

.00 

129.05 

.00 

.08 

13 

10.51 

.00 

.00 

186.06 

.00 

.18 

14 

11.49 

.00 

.00 

96.78 

.00 

.04 

15 

11.24 

.00 

.00 

102.21 

.00 

.16 

16 

14.5 

.00 

.00 

71.19 

.00 

.07 

17 

13.16 

.00 

.00 

141.80 

.00 

.09 

18 

13.99 

.00 

.00 

90.01 

.00 

.02 

19 

15.47 

.00 

.00 

81.01 

.00 

.12 

20 

13.49 

.00 

.00 

99.38 

.00 

.05 

21 

12.74 

.00 

.00 

136.40 

.00 

.13 

22 

15.25 

.00 

.00 

138.19 

.00 

.15 

23 

16.56 

.CO 

.00 

132.20 

.00 

.15 

24 

18.06 

.00 

.00 

100.18 

.00 

.25 

25 

17.70 

.00 

.00 

132.72 

.00 

.18 

26 

17.29 

.00 

.00 

142.58 

,00 

.17 

11 

16.94 

.00 

.00 

68.01 

.00 

.09 

28 

17.58 

.00 

.00 

69. 89 

.00 

.03 

29 

17.06 

.00 

.00 

41.89 

.00 

.19 

30 

16.31 

.00 

.00 

138.22 

.00 

.28 

31 

15.86 

.00 

.00 

149 i 30 

.00 

.01 

32 

18.16 

.00 

.00 

138.85 

.00 

.08 

33 

16.92 

.00 

.00 

148.49 

.00 

.11 

34 

15.21 

.00 

.00 

91.84 

.00 

.84 

35 

18.67 

.00 

.00 

75.41 

.00 

.52 

36 

18.48 

.00 

.00 

49.75 

.00 

.36 

37 

17.23 

.00 

.00 

89.32 

.00 

.16 

38 

18.14 

.00 

.00 

145.84 

.00 

1.36 

39 

18.12 

.00 

.00 

146.47 

.00 

.72 

40 

18.82 

.00 

.00 

160.19 

.00 

.51 

41 

19.22 

,00 

.00 

126.88 

.00 

.56 

42 

15.63 

.00 

.00 

162.42 

.00 

.90 

43 

18.47 

.00 

.00 

124.22 

.00 

1.12 

44 

19.44 

.00 

.00 

61.05 

.00 

.45 

45 

16.54 

.00 

.00 

135.41 

.00 

.84 

46 

18.17 

.00 

.00 

54.79 

.00 

.32 
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Table 5: Mode Shape Comparison (Concl) 

Case A * 100Z of Available Modes 
% Diff for Various Modal Coupling Tech* 


Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

47 

18.33 

.00 

.00 

142.35 

.00 

.30 

48 

18.25 

.00 

.00 

145.79 

.00 

2.22 

49 

19.16 

.00 

.00 

139.28 

.00 

.27 

50 

19*29 

.00 

.00 

65.342 

.00 

2.16 

51 

19.51 

.00 

.00 

62.72 

.00 

1.22 

52 

19.17 

.00 

.00 

52.42 

.00 

3.43 

53 

18.89 

.00 

.00 

105.50 

.00 

4.69 

54 

19.52 

.00 

.00 

38.21 

.00 

1.48 

55 

19.13 

.00 

.00 

103.39 

.00 

2.35 

56 

11.17 

.00 

.00 

161.11 

.00 

5.68 

57 

19.81 

.00 

.00 

11.29 

.00 

.17 

56 

19.81 

.00 

.00 

135,33 

.00 

.09 
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Table 6: Mode Shape Comparison 

Case B * 71X of Available Modes 



X 

Diff for 

Various Mddal 

Coupling Tech. 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWK 

RFIWOM 

*4 

14.84 

.00 

.00 

.00 

• 00 

.00 

5 

13.05 

.02 

.01 

.01 

.01 

.01 

6 

13.45 

.04 

.03 

.18 

.02 

.18 

7 

13.70 

.07 

.01 

.23 

.01 

.23 

8 

12,23 

.08 

,06 

.06 

.05 

.06 

9 

11.21 

.12 

.08 

.33 

.06 

.33 

10 

9.54 

. 16 

.18 

.16 

.15 

. 16 

11 

8.75 

.22 

.22 

.66 

.14 

.66 

12 

13.25 

.80 

.89 

.76 

.70 

.76 

13 

10.51 

.30 

.41 

1.78 

.27 

1.78 

14 

11.49 

.91 

1.21 

.88 

.57 

.88 

15 

11.24 

1.75 

2.05 

2.53 

1.67 

2.53 

16 

14.55 

.61 

.71 

.95 

.34 

.95 

17 

13.16 

1.02 

i.10 

1.02 

.70 

1.02 

18 

13.99 

1.52 

1.67 

1.63 

1.50 

1.63 

19 

15.47 

1.28 

1.03 

1.86 

.77 

1.86 

20 

13.49 

1.30 

1.03 

1.21 

.75 

1.21 

21 

12.74 

1.54 

2.02 

2.20 

1.23 

2.20 

22 

15.25 

1.93 

1.84 

1.98 

1.50 

1.98 

23 

16.56 

.38 

.62 

3.31 

.12 

3.31 

24 

18.06 

.32 

.81 

3.56 

.31 

3.56 

25 

17.70 

3.44 

3.63 

4.94 

3.43 

4.94 

26 

17.29 

9.71 

9.84 

8.91 

8.36 

8.91 

27 

16.94 

8.37 

8.24 

7.18 

7.04 

7.18 

28 

17.58 

7.16 

7.15 

6.09 

5.94 

6.09 

29 

17.06 

1.11 

2.22 

3.50 

.79 

3.50 

30 

16.31 

3.83 

5.48 

5.80 

3.52 

5.80 

31 

15.88 

3.90 

3.89 

4.95 

3.86 

4.95 

32 

18.16 

2.8) 

2.72 

4.45 

2.02 

4.45 

33 

16.92 

4.06 

3.29 

3.91 

2.92 

3.91 

34 

15.21 

6.36 

4.34 

25.34 

4.43 

25.34 

35 

18.67 

5.53 

5.12 

22.71 

.42 

22.71 

36 

18.48 

1.64 

1.62 

IS .07 

i . 79 

16.67 

37 

17.23 

4.65 

5.05 

12.42 

3.98 

12.42 

38 

18.14 

1.43 

2.71 

71.43 

6.38 

71.43 

39 

18.12 

2.45 

2.06 

189.04 

4.32 

189.04 

40 

18.82 

35.48 

19.60 

104.85 

60.42 

104.85 

41 

19.22 

39.52 

27.36 

128.81 

60.10 

J28.81 

42 

15.63 

26.78 

27.93 

126.75 

28.17 

126.75 

43 

18.47 

9.51 

9.90 

140.52 

38.29 

140.52 

44 

19.44 

14.36 

15.85 

144.04 

126.68 

144.04 

45 

16.54 

24.08 

24.65 

213.55 

142.78 

213.55 

46 

18.17 

3.25 

4.86 

150.22 

123.17 

150.22 
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Table 6: Mode Shape Comparison (Concl) 

Case B ■ 71% of Available Modes 



% D 

iff for 

Various Modal 

Coupling Tech. 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

47 

18.33 

4.54 

5.77 

250.51 

144.13 

250.51 

48 

18.25 

62.46 

70.50 

553.15 

77.72 

553.15 

49 

19.16 

57.97 

66.55 

965.00 

130.07 

985.90 

50 

19.29 

114.23 

116.94 

930.68 

128.01 

930.68 
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Table 7: Mode Shape Comparison 

Case C * 46% of Available Modes 
X Diff for Various Modal Coupling Tech. 


Mode Exact Freq 


No 

(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.02 

.02 

.03 

.02 

.03 

5 

13.95 

.18 

.04 

.03 

.03 

.03 

6 

13.45 

.47 

.45 

.49 

.31 

.49 

7 

13.70 

.31 

.10 

.55 

.13 

.55 

8 

12.23 

.54 

.23 

.18 

.16 

.18 

9 

11.21 

.1.59 

1.30 

2.00 

.90 

2.00 

10 

9.54 

1.10 

.74 

.41 

.36 

.41 

11 

8.75 

2.79 

3.01 

3.61 

1.79 

3.61 

12 

13.25 

2.50 

3.48 

2.42 

1.79 

2.42 

13 

10.51 

1.92 

2.35 

9.62 

2.83 

9.62 

14 

11.49 

4.84 

7.14 

5.19 

3.88 

5.19 

15 

11.24 

8.72 

9.37 

8.63 

5.92 

8.63 

16 

14.55 

3.33 

4.77 

3.71 

1.62 

3.71 

17 

13.16 

9.03 

90.46 

13.85 

5.90 

13.85 

18 

13.99 

11.17 

10.37 

13.28 

7.67 

13.28 

19 

15.47 

12.62 

11.13 

68.73 

9.12 

68.73 

20 

13.49 

13.23 

10.81 

81.34 

7.28 

81.34 

21 

12.74 

7.56 

16.01 

46.68 

10.40 

46.68 

22 

14.25 

16.60 

19.11 

54.62 

20.68 

54.62 

23 

16.56 

6.01 

5.57 

123.28 

97.95 

123.28 

24 

18.06 

4.87 

9.79 

129.15 

145.42 

129.15 

25 

17.70 

17.82 

16.45 

135.79 

155.65 

135.79 

26 

17.29 

68.71 

127.43 

150.09 

124.74 

150.09 

27 

16.94 

72.98 

144.11 

134.45 

120.45 

134.45 

28 

17.58 

103.74 

98.65 

152.44 

150.49 

152.44 

29 

17.06 

81.82 

139.33 

166.49 

128.91 

166.49 

30 

16.31 

110.88 

130.70 

145.30 

111.59 

145.30 

31 

15.86 

95.31 

153.24 

243.62 

141.77 

243.62 

32 

18.16 

101.15 

122.19 

695.80 

128.54 

695.89 
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Table 8: Mode Shape Comparison 

Case F * 27Z of Available Modes 



% 

Diff for 

Various Modal 

Coupling Tech. 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.04 

.05 

*06 

.03 

.06 

5 

13.95 

.24 

.07 

.26 

.06 

.26 

6 

13.45 

.71 

.89 

2.63 

.71 

2.63 

7 

13.70 

1.02 

.20 

8.69 

2.36 

6.69 

8 

r.23 

.74 

.39 

1.48 

.42 

1.48 

9 

11.21 

2.58 

2.75 

4.62 

1.96 

4.62 

10 

9.54 

2.84 

2.19 

1.40 

.98 

1.40 

11 

8.75 

10.05 

9.34 

38.19 

23.43 

38.19 

12 

13.25 

10.96 

11.05 

65.52 

28.35 

65.52 

13 

10.51 

7.03 

7.59 

128.49 

53.35 

128.49 

14 

11.49 

29.41 

33.79 

159.45 

43.11 

159.45 

15 

1) .24 

33.59 

66.95 

167.22 

88.27 

167.22 

16 

14.55 

30.24 

113.08 

161.12 

119.44 

161.12 

17 

13.16 

111.59 

138.72 

161.33 

140.08 

161.33 

18 

13.99 

123.87 

130.93 

148.55 

132.90 

148.55 

19 

15.47 

148.78 

127.29 

379.09 

127.01 

379.09 
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Table 9 

Comparisons of Maximum Absolute Values of Interface Loads 
Case A ■ 100% of Available Mo^es 


Load No. 

Exact Load libs) 

IMSL 

Percent 

I/F 

Difference - 
RFSWOM 

ABS. Max. 
RFIWM 

Loads 

RFIWOM 

1 

-481.999 

0 . 

-.00 

3419.56 

.00 

-.05 

3 

-202 . 1 38 

0 . 

-.00 

10154.67 

.00 

-.08 

5 

-498.819 

0 . 

.00 

3957.31 

-.00 

.11 

7 

474.713 

0 . 

-.00 

2208.71 

.00 

-.03 

9 

191.901 

0 . 

-.00 

22733.86 

.00 

-.02 

11 

486.870 

0 . 

.00 

17322.92 

.00 

.03 


IMSL * Inertial Coupling W/ Mass and Stiffness Loading 
I/F * Interface Method of Inertial Coupling 

RFSWGM * Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIVM * Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM - Residual Flexibility with Inertial Coupling, without Residual Mass 
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Table 10 

Comparisons of Maximum Absolute Values of Interface Loads 
Case B * 7 IX of Available Modes 


Percent Difference - ABS. Max. Loads 


Load No. 

Exact Load (lbs) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

i 

-481.999 

.10 

-.00 

2.08 

.10 

-.38 

3 

-202.138 

.26 

.27 

2.62 

.?■/' 

.26 

5 

-498.819 

.12 

.13 

-.13 

.14 

-.20 

7 

474.713 

.17 

.15 

-1.36 

.18 

.22 

9 

191.901 

.15 

.13 

.26 

.13 

-.24 

11 

486.870 

.19 

.19 

2.00 

.24 

-.10 


T MSL = Inertial Coupling W/ Mass end Stiffness Loading 
i/F * Interface Method of Inertial Coupling 

RFSWOM * Residual Flexibility with Stiffness Coupling- without Residual Mass 
RFIWM * Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM * Residual Flexibility with Inertial Coupling, without Residual Mass 
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Table 11 

Comparisons of Maximum Absolute Values of Interface Loads 
Case C * 46X of Available Modes 

Fercent Difference - ABS. Max. Loads 


Load No. 

Exact Load (lbs) 

IMSL 

l/F 

RFSWOM 

RFIWM 

RFIWOM 

i 

-481.999 

-1.32 

-1.29 

-4.07 

-1 .09 

.29 

*> 

jf 

-202. 138 

2.34 

2.33 

3. 10 

2.43 

).I9 

5 

-498.819 

.23 

.39 

-7.33 

.43 

. 1 7 

7 

474.713 

-1.48 

-I .43 

2.93 

-1.36 

-.34 

9 

191.901 

.00 

12 

-23.44 

-.07 

-1.48 

11 

486.870 

-.83 

-.77 

-3.83 

-.78 

-.7 1 


IMSL - Inertial Coupling W / Mass and Stiffness Loading 
I/F ® Interface Method of -nertial Coupling 

RFSWOM * Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIVJM * Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM * Residual Flexibility with Inertial Coupling, without Residual Mass 
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Table 12 

Comparisons of Maximum Absolute Values of Interface Loads 
Case F ” 27Z of Available Modes 

Percent Difference - ABS. Max. Loads 


Load No. 

Exact Load (lbs) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

i 

-48 1 . 999 

-2.45 

-2.88 

6.38 

-6.66 

-2.60 

3 

-202.138 

3.04 

-2.85 

5.68 

-10.99 

.27 

5 

-498.819 

-2.24 

-2.61 

8.02 

-4.49 

-1.85 

7 

474.713 

-3.82 

-4.44 

10.81 

-8.79 

-10.31 

9 

191.901 

-2.79 

-.54 

21.52 

- 7.54 

-14.78 

11 

486.870 

-.08 

-.48 

1 7 .06 

-4.28 

-7.60 


IMSL = Inertial Coupling W / Mass and Stiffness Loading 
l/F ~ interface Method of Inert ini Coupling 

RFSWOM = Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIWM = Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM = Residual Flexibility with Inertial Coupling, without Residual Mass 
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Another comparison study of modal synthesis methods was conducted by 
W. Benfield , C. Bodley and G. Morosow. Here are a few extracts 
from that publication. 

Some of the earliest formal documentation on modal synthesis is 
accredited to Hunn and Gladwell. Although their methods yield good 
accuracy, they are restricted to structures with statically 
determinate interfaces and have been somewhat overshadowed by more 
recent developments. In this paper, only methods applicable to 
redundant interface connections are compared. 

Hurty developed the first modal synthesis method capable of analyzing 
structures with redundant interface connections. In this method, the 
component vibration modes are determined with all interface 
coordinates fixed, and are thus callled fixed-constraint modes, and 
fixed-constraint, natural vibration modes are used as the generalized 
coordinates to determine the modes of the system. 

Bajan and Feng and Craig and Bampton introduced the concept of using 
only constraint modes, thus eliminating the need to distinguish 
rigid-body modes and redundant constraint modes. The 
fixed-constraint, natural-vibration modes of each component are 
determined in the same manner as in Hurty* s method. The constraint 
modes and fixed-constraint, natural-vibration modes are then used as 
the generalized coordinates to determine the system modes. 

Goldman developed a method of vibration analysis using dynamic 
participating. In his method, the component vibration modes are 
free-free, rather than fixed-constraint modes. Rigid body modes and 
selected component vibration modes are used to solve the 
interconnection problem. This eliminates the need for constraint 
modes. The component free-free vibration modes are then used as 
generalized coordinates to determine the modes of the system. 

Hou has described a method for medal synthesis that also uses 
free-free component vibration modes. Constraint modes are not used. 
Selected component modes are used to connect the components 
together. These modes then become dependent coordinates thus 
reducing the number of independent, generalized coordinates by an 
amount equal to the number of system interface coordinates. 

Benfield and Hruda presented a method of component mode substitution 
using either free-free or constrained component modes called branch 
modes. Constraint modes are used for constrained components, but 
these are reduced to dependent coordinates rather than carried in the 
solution as independent coordinates. Only the component vibration 
modes are used as generalized coordinates to determine the system 
modes. In addition they also presented the concept of using 
interface-loaded component modes to improve the accuracy of the 
solut ion. 
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In general, the methods that used constrained-component vibration 
modes produced more accurate results than those that used free-free 
vibration modes. Furthermore, methods that used component mass and 
stiffness matrices were more accurate than those that did not. 

The features of the methods can be summarized as follows. 

Method 1 - Benfield and Hruda, Free-Free Component Modes 

- Less accurate 

- Uses free-free component modes 

- Easy to use, since no constraint modes or matrix inversions are 

used. 

Requires mass and stiffness matrices for the components. 

Method 2 - Benfield and Hruda, Free-Free Component Modes, with 
Interface Loading 

- Good accuracy 

- Uses free-free, interface-loaded component modes 

- Requires constraint modes or reducing transformations for each 
component • 

- Requires mass and stiffness matrices for the conponents. 

Requires interface loading of component modes. As a result, 
other attached components affect the component modes. 

Method 3 - Benfield and Hruda, Constrained Component Modes 

Good accuracy 

- Uses constrained-component branch modes. 

More difficult to set up for many components 

- Requires constraint modes for each component 
Requires mass and stiffness matrices for the components 
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Method 4 - Benfield and Hrude, Constrained Component Modes with 
Interface Loading 


One of the moat accurate 

Uses constrained-component branch modes 

More difficult to set up for many components 

Requires constraint inodes for each component 

Requires mass and stiffness matrices for the components 

Requires interface loading of component modes* As a result, 
other attached components affect the component modes* 

Method 5 - Hurty 

" Ona of the most accurate 

Uses fixed interface component modes 

Requires mass and stifness matrices for the components 

System rigid-body and redundant interface generaiited 
coordinates must be retained as additional degrees of freedom in 
the modal solution* For structures having a large number of 
interface coordinates, the number of component modes used may be 
limited. 

Method 6 - Bajan and Feng, and Craig and Bampton 
One of the most accurate 
Uses fixed-interface component modes 

Requires constraint modes. This method is more convenient to 
use since it is not necessary to identify statically determinate 
and redundant interface coordinates. 

- Requires mass and stiffness matrices for the components. 

System interface generalised coordinates must be retained as 
additional degrees of freedom in the modal solution* For 
structures having a large number of interface connect iors , the 
number of component modes used may be limited. 



Method 7 - Goldman 


One of the least accurate. 

Uses free-free component modes. 

Modal selection is important. A matrix inversion is required, 
and the matrix must have the same size as the number of 
interface coordinates. Since an ill-conditioned inversion may 
result from improperly selecting the component modes, it is 
desirable to have as few interface coordinates as possible. 

The modal stiffness matrix is nonsymmetric 

Does not require mass and stiffness matrices for the 
components. Thus, system modes may be synthesized using only 
the component modes. 

The first n modes in the system eigenvalue-eigenvector solution 
are not usable. (Here n is equal to the number of system 
interface connections). 

Method 8 - Hou 


One of the least accurate 
Uses free-free component modes 

Modal selection is important. A nonsymmetric matrix inversion 
is required, and the matrix must have the same size as the 
number of interface connections. An ill-conditioned inversion 
may result from improperly selecting component modes. Thus, it 
is desirable to have as few interface coordinates as possible. 

Does not require mass and stiffness matrices for the 
components. Therefore, system modes may be synthesized using 
only the component modes. 

The number of component modal coordinates used in the system 
modal solution is reduced by n (where n is equal to the number 
of system interface coordinates). 



B. SHORT-CUT TECHNIQUES 
1. Introduction 


Jn Section A we discussed sever*! techniques to obtain the response 
of a booster/payload system. All these techniques necessitate the 
solution of the coupled booster/payload equations, i.e., they are 
“full-scale 11 methods. As we pointed out before, this solution is quite 
expensive, especially if it has to be repeated several times, for example 
during a design effort. Although mass and stiffness changes during a 
design effort are often small, current practices used to design ^ayload 
structures require a new "full-scale" solution every time such small 
changes in the payload are made. A similar situation exists in the case 
of payloads that are designed for multiple flights with moderate 
configuration changes. 

A need exists for the development of "short-cut" methods. The term 
"short-cut" method implies that the method should be able to evaluate 
small changes in the payload in a relatively short time. First, a 
short-cut method should avoid the direct solution of the coupled 
equations of the booster/payload system. Secondly, it should avoid as 
much as possible the interfacing between different organizations. This 
means that one should strive towards as much independence for the payload 
design organization as possible. 

The objective then of Section B is to present several of the most 
promising of these short-cut methods. The first of these methods will be 
discussed in the next section, 

2. The Perturbation Technique 

In this section we shall discuss a short-cut method which is based on 
a well known perturbation technique. We shall first discuss the 
perturbation technique in general terms and then apply it to the 
particular problem of a booster/payload system. 

Let us consider a set of equations of motion of a certain structure, 

[» 0 ] i * I + [k 0 ] } x ! - | F [ now 

where [Mq] and [Kq] are the mass and stiffness matrix of the 
structure, respectively. The vectors | x} and ^f} are the generalized 
discrete displacement and force vectors. The eigenvalue problem 
associated with Equation (104) can be written as 

(-“o ["ol + [»„] ) |»ol ' M (105) 
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The solution of this eigenvalue problem yields a modal matrix [^q] and 
a diagonal eigenvalue matrix f *i satisfying 


[*o] ["oHVl 


t‘3 • w'tgtg 



(106) 


Next, let us assume that the elements of [Mq] and [Kq] undergo 
small changes, so that the new system of equations can be written as 
follows 


[M]jx| ♦ [out - M 

with 

[M] - [Mj t C[MJ , [k] - [k q ] + r[ Kj ] 

where £ is a small parameter such that 

E [Mj * [M] - [»<,] . E 0,] - M - [? 0 ] 


(107) 


(108) 


(109) 


Note that the matrix differences on the right-hand sides of Equations 
(109) are small, so that it is easy to determine a small £ so that [MjJ 
and [Kj] are of the same order of magnitude as [M] , [Mq] and [K], 

[Kol. 

The objective of the perturbation technique is to obtain a solution 
for the new eigenvalue problem associated with Equation (107). 

[m] + [K] ^ \ *\ = {of (1 

without actually solving Equation (110). To this end, let us write 


w - 

W + 

e \ x l\ 

+ e 2 { x 2 } + 

dm 

1*1 . 

{» 0 f + 

C W + 

e 2 ) <f 2 } + ... 

(112) 

(i) = 

“o + 

e Ml 

2 

+ e (x >2 + 

(113) 


ORIGIN . r u- 
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This can be done because of the small changes in [Mq] and [Kq] as 
expressed in Equation (108). Also, 


M " + e + e 2 { q 2 } + 


(114) 


where { qoi an< * { q} are the normal coordinate vectors of the unperturbed 
and perturbed system, respectively, i.e. 


I* I * O ] 1 it 
K I ’ [* 0 ]ti 0 > 


with 


(115) 


(116) 


0 .1 T [«][♦] -t I] , [*] T [K][*] t [K][*] * C» 2 3 (117) 

First, let us substitute Equation (116) into Equation (107) and 
premultiply by [<b] T yielding 

|q| + *\ = C^] T | F} (118) 

Substituting Equations (112), (113) and (114) into Equation (118) and 
equating coefficients of like powers oft* , we can write, 

T 

+ K} = [♦„] i F f (H9) 


W + [%] Kl = [*i] T \?\ - 2 OolOJ M O20) 

It is now possible to solve Equations (119, 120, etc.) in sequence. 

The f:*cat of Equations (119) represents the unperturbed equation of 
’TtOt ion, i.e., the modal form of Equation (104). This solution is 
avai able or can be determined. 

Once the vector {q 0 } is determined one can solve Equation (120) if 
[4^] and [<^] are known. The determination of these matrices is the 
subject of the next paragraph. 


ORtcyti 

OF POOR 


1. 


■*» I* 

■ I f 
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First, it is always possible to write as a linear combination 

of the eigenvectors [4 >q] , 

OJ = M f a ] 

( 121 ) 


+ 

+ 


because ^q] is a complete set of vectors (i.e., they can be used as a 
basis for a vector space). Note that [cx ] represents the coefficient 
matrix of ] in the linear combination and must still be determined. 
To this end, let us introduce Equations (108), (112), (113) and (121) 
into Equation (117), and only keep terms inlands', 


c[ *<j w) (ftj * E [")]) ([♦„] + 

£ ftj 0]) ( t K o] + £ [K,] ) ( W + E 


e [♦<)] [“])- [l] 

M [«] ) , , 

([%)]+ *,]) 


(122) 


(123) 


j' | 

Equating coefficients of like powers in S' and £ we obtain from 
Equat ions (122-123) 


M 

2 K] [ w i] - 


+ £0 T - - [^o] T Oil M 

t“o]M + - KH 


(124) 


(125) 


Equations (124-125) can now be solved for and [U^]. This 

enables use to solve Equation (120) for {q^J , so that from Equation 
(114) we obtain the first order approximation 

{<?} = |q 0 } + e {qj 0.26) 


and from Equation (116) we obtain { x} where we use 

[♦] = M + r OJ 


(127) 


This perturbation technique can now be applied to a booster/payload 
situation. The assumption is that only changes of order £ are made in 
the payload, i.e. , 


W = + f - [MpJ ’ W = + 1 [Kpil 


(128) 


r~ 
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where [Mpg] and [Kpg] are the mass and stiffness matrices of the 
unperturbed payload Pq. Let us write the mass and stiffness matrices 
in Equation (25) for the perturbed payload P. 


_m b + t ?Vp| 


k b + t ?Vp ! 0 

I P M P T P j I p M p I p 

» 

° 1 T pVp _ 


with 


[ t p] 



Using Equation (128) we can write [Sp] as 


( 129 ) 


(130) 



(131) 


and, keeping terms in £ only, we obtain 


Pp] 




(132) 

[C]) 


or 

Pp] ■ [ S P 0 ] + E OpJ 


(133) 


with 



(134) 



(135) 


It is now possible to write [Tp] as 


[Tp] - 

Ppo] + 6 [ T Pl] 

(136) 

with 



[ t p,j ■ 

[° H . h PI ] - [° S »1 

L» o J 1 L') i J 

(137) 


C’.'-'G,?:':. . _ 
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Finally we can write the mass and stiffness matrices of the coupled 
booster/payload system as follows 




‘ t«ciJ ' M M + 



(138) 


where C indicates Coupled. Matrices [M^] and [K^] are given by 
expressions (129). Taking into account Equations (128) and (136) we can 
write 
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and 
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where [Tpg] and [Tp]J are given by Equation (137) • Equation (138) is 
now equivalent to Equation (108) and the perturbation technique can be 
applied. 


Note that theoretically one can also obtain the higher-order 
perturbations , t , etc. But for all pratical purposes only £ 
perturbations are included. The question then is, how important , £ 

, etc. terms are. It is evident that Equation (111) is only va^id as 
long as the asymptotic expansion in £ does not break down, i.e., as long 
ascl^it^^ j etc. There are indeed cases where such an 

asymptotic expansion breaks down. It is then necessary to introduce 
other perturbation techniques (e,g. Lighthill, multiple scales, etc). In 
this connection it is important to recognize the fact that small changes 
in the mass and stiffness of the payload produce small changes in the 
eigenvalues but not nece ssar ily in the nodes. 
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3. The Base Drive Technique 


ORifiiiVAL r,u __ -.j 

OF POOR QUALITY 


In this section we shall discuss the Base Drive Technique as 
developed by k . Devers, et.al* Let us first recall Equation (91) in 
partitioned form, 
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where we solved for j in Equation (144) 

The payload designer is primarily interested in predicting the 
response of the payload which is given by 

f 4 ‘ 


(144) 


with 
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(145) 


where {a£} and {xl^jmust be computed from Equations (142-144) and 
from Equation (90;. 

'* R) 

The idea of a base drive short -cut method is to approximate {*£ j in 
Equation (144) in uuch a way as to avoid the solution of the complete set 
(142-144). To evaluate a particular short-cut method, the approximation 
of {x^j must be compared to the exact value given by Equation (144). 

A significant simplification of Equation (144) occurs when the 
interface is statically determinate, i.e., when 


-p 1 

_^N_ 


-P 

X N 


1 * 


«K 

liR 

X I , 


x r 

' *1 


(146) 




Kvb] • [ t Jvp] ■ f°] 


(147) 


This 
becomes , 


eliminates the dependence of 
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A first possibility is to assume that the presence of the payload has 
no effect on the response of the booster, i*e. $ = {0j* This 

approach is ca:ied the Direct Base Drive Technique or Open Lo~p Base 
Drive. Indeed if then Equation (143) becomes 

DpVp] + [^Vp] r4> ■ - ['pVpj lx f 1 aw) 

which means that the payload is "directly" driven at Its base (i.e. its 
interface with the booster B) by the force on the right hand side of 
Equation (149). Theoretically, the Direct Base Drive Technique assumes 
there is no coupling at all t between the booster and the payload. 
Practically, it means thatfj^ J^iX^for that the feedback of the payload 
is negligible. The conditions under which such an approximation is valid 
is still an unanswered question. The development of a criterion of 
validity of the use of the Direct Base Drive Technique should be 
considered. 


4. The Impedance Tec h nique 

In this section we shall discus yet another approach to the solution 
of the equations of motion of the booster/payload system. The Impedance 
Technique as developed by K. Payne is basical a f 11-scale method in 
the sense that it does not make any assumptions concerning the size of 
the payload nor the extent of the changes made in the payload. However, 
the method does avoid a full-seal c solution of the coupled 
booster /payload equations of motion and is particularly suited to deal 
with small changes in the pay load. The Impedance Technique is 
essentially a Base Drive method (see Section 3). It differs from tne 
approach in Section 3. in the manner in which the interface 
accelerations £ xfj are computed. Indeed, the interface 
accelerations will be computed in the frequency domain instead of the 
discrete time domain thereby essentially converting c set of differential 
equations into a set of algebraic equations. 


Let us now derive tie necessary equations. 
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First, recall equation 



(150) 


and write the top and bottom partitions separately, 
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where we invoke' 4 Equation (2) and assumed the external payload fore*** 


to be absent 
Equation (87) • 


(l F P f - {of)* we intro^jce an equation similar to 


E p 

i x B l - I*bI +)x b( 


(153) 


v lere the F- vector represents the booster response due to the externally 
applied force vector {Fj^and therefore satisfies, 

[“bJki +[<]m - \ (i54) 

and the R-vector represents the response of the booster due to the 
feedback of the payload and satisfies 

["»]<*»( R + [ k b]) x b | R -{-«?| < 155 > 

Applying transformation (32) to Equation (154) yields. 




(156) 


We now consider Equation (153) and recall Equation (4), 



Keeping the bottom partition in Equation (157) yields 


(157) 


jx® (t)j - jXj (t) | + jx® (t)j (158) 

where we now included the dependence of the vectors on the discrete 
time t. Taking the Laplace Transform (with zero initial conditions) of 
both sides of Equation (158) introducing the transformation s sjil (with s * 
Laplace variable * the ith input frequency; and j - r=9 we can 
write 
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which represents Equation (158) in the frequency domain. Taking the 
second time derivative of Equation (159) yields: 


i**n 

1*1 ««!> 
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Ihe basic idea of the Impedance Technique is to calculate the 
interface accelerations in the frequency domain and then 

transform them back to the discrete time domain, Ihe two terms on the 
righthand side of Equation (160) will, be replaced by algebraic matrix 
expressions so that the calculation of does not involve the 

solution of a set of differential equations. Let us start with the first 
term on the right-hand side of Equation (160), To this end let us 
convert Equation (156) to the frequency domain 



or, taking the second time derivative, 
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Let us write Equation (32) 



with 
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Then the bottom partition of Equation (163) 
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Premultiply Equation (162) by 



and invoking Equation (165) yields 



(J0 t ) 
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Equation (167) then, yields the first term on the right-hand side of 
Equation (160). The matrix A(jjl,) in Equation (168) is the transfer 
admittance from the points of application of the external forces {Fg} 
to the interface accelerations. 

Similarly, Equation (155) can be transformed into 
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where this time 
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is the matrix of coefficients for the point admittance for the booster at 
the interface. Equation (169) yields the second term onJthe right-hand 
side of equation (160). However, the reaction vector £R£(j A*) is not 
known a priori. To determine this vector let us consider Equation (152) 
which represents the payload equations of motion, and write it as 
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where we used Equations (5) and (7). Introducing the modal 
transformation (39) and taking into account the properties (40) we obtain 
from Equation (171), 
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with 

H] ■ M M 
M ■ W W * W 

r pi 

and L given by Equation (28). Note that I I zer ‘'hen 
interface is statically determinate. L J 
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The top 


and bottom partitions of Equation (172) can be written as, 
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ok: 

of kk; v- 

He shell now usuw Chet Che interface ie statically 
deterainate/L<f Ja{o])and calculate an expression for {r£} frosi Equations 
(175-176). First, we transfora Equation (175) to the frequency domain. 
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and fron Equation (176) ve obtain. 



Substituting Equation (177) into Equation (178) yields .* 
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from which we obtain the following expression for 
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is the impedance matrix of the payload at the payload/booster 
interface. Finally, we substitute Equation (180) into Equation (169), 
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Combining Equations 

(160), (167) and (182) yields 
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where we also used Equation (1). 

The coefficient matrix of £^ C 2 r(/^/)j , in Equation (183) represents the 
coupled impedance of the booster/payload system, and the right hand side 
represents a pseudo-generalized force. The interface acceleration can 
now be computed from Equation (183) with relative ease. 
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If we now consider a new payload on the same booster and with the 
same force {Fg}, the right-hand side of Equation (183) does not change 
so that , 

m\ 
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provided the interface does not change. 
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The interface accelerations in the discrete time domain can now 

be derived from equation (183) or Equation (184), using the inverse 
Laplace Transform. The payload response then follows from Equation (173). 


The approximation involved in the Impedance Technique is imbedded in 
the transformation to and from the frequency domain. If these 
transformations were exact, the method of determining {Xj-ftfy would be 
exact. Therefore, a detailed investigation of these transformations is 
required. There are also some problems pertaining to the modal damping 
when working in the frequency domain. 


Although Equations (183) and (184) were derived for an undamped 
statically determinate system, it is clear that damping and statically 
indeterminate interfaces can be included. For an indeterminate system it 
becomes necessary to keep track not only of the interface accelerations 
but also of the velocities and displacements at the interface. The use 
of the Fast Fourier Transform in obtaining the spectral data to be used 
in Equations (183-184) presents some problems. In general, however, 
enough correlation with the exact time domain solution is apparent to 
warrant further investigation into possible improvements. 
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5. The Generalised Shock Spectrum Technique 


The approach presented in this section is a generalised version of 
the shock spectrua technique as developed by Bamford. 

In order to explain the basic ideas underlying this technique let us 
recall Equation (43) 


'•♦KWi 1 1 1 "b 


♦f'XVB ! 
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(185) 

and let us assume we retained M modes for ^he booster (i.e.O^is an MxM 
matrix) and N modes for the payload (»e- ] is an NxN matrix) 

The basic idea of the shock spectrum technique is to determine load 
maxima without having to solve Equation (185). To reach this goal, a new 
model both for the coupled system and the forcing function ^Fg^ is 
introduced. 

First, the (N + M) modally coupled Equations (185) are replaced by 
(N x M) sets of two simultaneous equations each of which represents the 
coupling of one payload mode with one booster mode, as follows: 
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where we assumed that the interface is statically 
determinate ( i C- C“ r p K p" 7 p3 - CO ] ) 

Secondly, a bound qgp on each of the (N x M) modal responses of the 
payload is established. This is done by introducing a new model for the 
forcing function in Equation (186). The rather complicated forcing 
function is replaced by a much simpler function (e.g», an impulsive 
force) which produces the same maximum response peak as the original 
force. Therefore an analytical solution for both the response and 
maximum response of Equation (186) j.s possible (after some additional 
simplifications). Finally, a bound qp on the total payload response is 
constructed by summing over all the individual modal bounds qgp (over 
absolute values or in a root -sum-square sense that allows for phase 
weighting). Payload member loads are obtained by adding the 
contributions of all payload modal loads. 

As stated above, the forcing function in Equation (186) is replaced 
by a modal delta function of a certain magnitude Fg. This magnitude 
Fg is evaluated from an already existing transient analysis of the 
booster with or without a dummy payload. 

The main objections that can be raised against the Shock Spectrum 
Approach are: 

1) No critical evaluation is available on the validity of replacing 
model (185) by model (186). What effect does this replacement 
have on the load bounds? This change of model could not only 
result in a toe conservative design but also in an 
unconservative one. Model (186) not only ignores the coupling 
between the B-modes due to rigid body feedback of P , but more 
importantly it ignores the effects that the coupling of one 
B-mode with one P-mode has on all the other P-modes. 

2) The manner in which Fg is calculated again leaves the question 
of whether or not the envelope values are conservative or not 
and by how much. 

3) The method appears rather complicated and is not simple to use. 
This can lead to misinterpretation and confusion when the method 
is applied. More rigor in the mathematical formulation is 
desirable. 

The method can be of definite use in the first stages of a design 
effort. It leaves the payload designer the freedom of operating 
independently from the booster organization. This results in a quick 
turn around time. The method could then in the final stages be 
complimented by a more rigorous approach. 
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CHAPTER II- A PAYLOAD INTEGRATION TECHNIQUE - 
FULL-SCALE VERSION 


1- Introduction 


Present analytical techniques by which design loads are predicted are 
very costly and time consuming. The calendar turnaround time of a given 
cycle usually is lengthened when the payload design organization, the 
booster organization and the payload integration organization are 
different companies. Indeed, a fair amount of coordination is necessary 
to make the transfer of information between those three organizations 
optimal. Unfortunately, this coordination is very difficult to 
establish, resulting in considerable time delays. Moreover, these costs 
and delays repeat themselves for every load cycle (i.e., every time a 
change is made in the booster or payload). A typical example is the 
development of the Viking Orbiter System as mentioned in Chapter I. 

The ever increasing number of modal coordinates necessary to model 
today's aerospace structures not only increases the cost of a load cycle, 
but also imposes greater demands on the analyst to keep the models within 
range of current computer capabilities. 

The objective of the present chapter is to develop a "full-scale" 
payload integration method which reduces the cost of a load cycle and 
will be capable of handling ve*y large systems. This new approach is a 
"full-scale" method in the sense that it actually solves the coupled 
booster/payload system equations and does not involve any additional 
approximations or assumptions as compared to the standard transient 
analysis. 

The configuration of multiple payloads connected to the booster 
through separate interfaces is typical for most Shuttle missions. The 
fact that many of these payloads are not directly coupled allows for a 
significant simplification of the booster/payloads system equations. 
Superfluous interface degrees of freedom on the booster side can be 
accommodated within the formulation. Superfluous interface degrees of 
freedom are those interface degrees of freedom which are included in the 
launch vehicle model but are not connected to any payload. The 
superfluous interface degrees of freedom arise because the booster 
organization cannot afford to reconstruct a booster model every time the 
interface with payload(s) changes. 

A numerical integration scheme used to obtain the booster/payload(s) 
system response is defined. The standard approach is to obtain the 
so-called "modal modes" i.e. the coupled system modes in order to 
decouple the system equations. The present approach avoids the solution 
of such a system eigenvalue problem. A Newmark-Chan-Beta numerical 
integration scheme is used to directly determine the system response. 


69 



This technique takes advantage of the peculiar structure of the equations 
of motion for the system* All zero entries of the system mass and 
stiffness matrices can be eliminated. In fact, a comparison feature can 
be implemented so that elements close to zero can be omitted, reducing 
the cost of the integration routi .e even more. The comparison feature 
makes the full-scale method a so-called short-cut method. A full-scale 
method can be defined as a method which yields "exact 11 results compared 
to the standard transient analysis technique, whereas a short-cut method 
introduces an approximation or assumption which leads to a more cost- 
effective but less accurate solution. 

A Fortran computer program has been written and implemented on the 
CDC Cyber 172. The final remarks section discusses advantages and 
disadvantages of the proposed approach. 

2. Derivation of the Basic Equations 

The objective of this section is the derivation of the equations of 
motion of the launch vehicle/payload(s) system. Figure 3 shows the free 
body diagrams of the booster B and the payload PI and P2. The booster 
and the payloads are connected to each other through the interface. We 
also consider superfluous interface coordinates. Indeed, often the 
booster interface contains more degrees of freedom than is necessary to 
attach the payloads. The reason for this is that the booster 
organization provides a set of interface restrained booster modes which 
can be used to accommodate many different payload configurations. It 
would be prohibitive to recalculate a set of cantilevered launch vehicle 
modes every time the interface with the payloads changes. 

From the free body diagrams in Figure 3 we can easily write the 
uncoupled equations of motion for the booster B and the payloads PI and 
P2 (This is similar to equation (3)): 
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Figure 3 Free-Body Diagrams of Booster B 
and Payloads PI and P2 






Several well known procedures to couple the above three equations 
were discussed in part A of Chapter I. Again, the objective is to 
eliminate the redundant interface displacements and, in the process also, 
the unknown reaction vectors. One such technique uses so-called 
"cantilevered" booster and payload displacements. This technique was 
discussed in section A-3f of Chapter 1. To set the stage, let us note 
that we can partition the vectors £xgf , |xpif » | x p2f » i » 
etc. as follows, 



Similar partitions can be written for the velocities and the 
accelerations. In the above equations we assumed that the external 
forces only act at the non-interface booster degrees of freedom. This 
assumption is only made for convenience. Also It should be noted that 
this development is not limited to two payloads. 
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g - pi p2 

Next, it is easy to show that i x ^f • and ? x / V J can be 

written as follows: 
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These equations are similar to Equation (81). 
Using Equation (189) ana the fact that 
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we can form the following transformation, 
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where , 
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NOTF: 

IF « IF1 + IF2 + IFS 
B - KB + IF 
Pi - NP1 + IF1 
P2 * NP2 + IF2 

also, 0^, the zero matrices 0^, 0 0^, 0^, 0^, 0^ 0^, 

°8» 0 9> 0 1Q , 0 11 *^ 12 » °i 3 » and 0 i 4 have dimensions BxNPl, 

BxNP2, NPlxIF2, NPlxIFS, PlxNB, PlxNP2, IFlxIF2, IFlxIt’S, NP2xT.Fl, 

NP2xIFS, P2xNB, P2xNPl, IF2xIFl and IF2xIF2 respectively. 
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Taking into account that the reactions between interfaces are equal 
but opposite we can substi :e transformation (191) into Equation (187) 
and then premultiply by the transpose of the transformation matrix. This 
will eliminate the redundant interface displacement vectors as well as 
the reaction vectors. The resulting equations are the couple discrete 
equations of motion and cen be written as: 
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where , 
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and again, it can be 3hown that [I'gKgTg], f^pi K pi T p)J snd 
[ I p 2 K F 2 T p 2 J alw '»y s vanish and iT^Tg], [Tj 1 K p 1 T pl l and 
^^P2 K P2 T P2^ vap ^ stl when the corresponding interface is determinate. 


Next, we c .reduce the interface restrained modes, i.e. 
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These mo^ s :an be truncated ccording to a predetermined cut-off 
frequency, hereby reducing the si*re cf Equation (193). Substituting 
Equations ; 5; into Equations (193) we obtain: 
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Note that we introduced damping and that we kept the interface 
degrees of freedom in discrete forsK This is recommended when DENSE 
Logic is used in programing the procedure. If PARTITION Logic is used, 
we actuall, introduce the interface modes l^jli 


f£‘K M ii]#>i] - hj 

#>>^1 - tco*i 


Although the booster and payload modes can be truncated, it is imperative 
not to truncate any interface definition. This is of special 
significance if one introduces the interface modes. The motivation for 
this will be explained in the load section. 


Before proceeding with the next section, we wish to introduce the 
following notations: 
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So Chat Equation (196) can be written as: 

iM] {q} + [C] {q£ + [K] {q} = \ f } 
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3. The Numerical Integration Scheme 

The objective of this section is to choose a suitable integration 
schema to directly integrate Equation (201). This is in contrast with 
the usual approach, where a new eigenvalue problem with [M] and [K] as 
mass and stiffness matrix is solved. 

It is important to recognize that system (201) can have a frequency 
content much higher than the cut-off frequency of the forcing term. 
Usually, those higher frequencies will not produce significant 
responses. Therefore, a suitable numerical technique should be capable 
of using a stepsize h, which reflects only the highest frequency of 
interest but at the same time remains numerically stable. For example, a 
Runge-Kutta routine would not be suitable because it requires a time-step 
consistent with the highest frequency in the system, even though these 
high frequencies tiay not be of interest to the analyst. Although there 
are techniques to obtain a good estimate of the highest system frequency, 
using such a frequency to determine the stepsize would unnecessarily 
increased the cost of the response routine. 

A method that satisfies above constraints is given by the 
Newmark-Chan-Bet a integration method: 


{q> i+1 = (qh + (1-Y)h{q) i + Y h{qh +1 


( 202 ) 


«>ifl * «>i + Mli + (0.5-wn 2 «) 1 'fh 2 {q} 1+l 


+ [CJ (qj 1+ l + («{<,} 1+1 - 


where the method is unconditionally stable if (23^+1)2/16. 

Artificial positive damping is introduced when 0.5, and artificial 
negative damping if < 0.5. Good values for our purposes are 0.5 
and ^9* 0.25. Theoretically, the time step h can then be given any value 
while the scheme remains stable. In fact, for very large values of h, 
the scheme generates the static solution of Equation (201). \l 3 c, the 

scheme will damp out the highest (and least important) mc^s while 
preserving the lower ones. In addition, as we shall are shortly, the 
Newmark-Chan-Bet a scheme is capable of taking advantage of the peculiar 
structure of the present equations of motion. Indeed, let us sub.^itute 
Equations (20.?) and (203) into (205) and obtain 


(203) 


{f} 


i+1 


(204) 
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[M+*hC+/3h 2 K]{q$ i+1 «{ F } i+ l - 

- [ (l-T)hC + (ij -^)h 2 K] {q} t 


[C + hK] {4} t 

- lK]{q} i 
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or, using Equation (199), we can write: 
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are diagonal matrices, except D 2 , D 5 and Dg which are IFxIF 
matrices. Consequently, the evaluation of fg£, fj£, and fp£, are 
very cost-effective. Also, note that Dj (i*l,...,9) are one-time 
calculat ions . 


■ t \ c *-• 

(205> 

(206) 

(207) 

(208) 
(209) 
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Normally, the soluticu for Jq} in Equation (206) requires a 
triangular decomposition and must be repeated for every h* However, in 
this case the unique form of the coefficient matrix of £q J £+1 makes it 
possible to avoid such decompositions* Moreover, it is possible to 
completely take advantage of the diagonal and zero partitions appearing 
in that coefficient matrix* 

First, let us premultiply Equation (206) by {-BDj“l j I | -PDjl ) 


which yields the following expression; 

(*Ii*ll -AxWui ♦ fli ♦ A 3 f pi } (211) 

with 

Ai - ( D 2 + A,B T + A 3 P T ] _1 (212) 

A 2 - - BD " 1 , A 3 - - PD ' 1 (213) 

Furthermore, from Equation (206), we easily obtain, 

*Ni*l * D'^Bi ♦ Apf ? i+1 ( 214 ) 

^Si +1 " ° 3 lf Pi + ( 215 ) 


Equations (211), (214) and (2lS) represent the final set of 
recurrence relations replacing Equation (206)* Note that Equation (212) 
represents the inversion of an IFxlF matrix where IF is relatively small 
in many applicat ions. This effectively removes the problem of triangular 
decomposition of a large matrix* Also, note that the cost of the 
algorithm primarily comes from multiplications involving matrices h\, 

Ao and A 3 * Note however, that their dimensions are IFxlF, IFxNB and 
IFxNP respectively* In addition, the routine requires much less core 
memory which allows for the so* .tion of much larger problems* 

4* The Load Transformation 


The purpose of this section is to briefly review the "accelerat ion 1 ' 
approach to calculating loads and at the same time point out some 
possible savings. An elementary member load transformation can be 
written as: 

{t-j} • [kVlj (*pji ( 216 ) 

Note that such a member could he part of any payload. For example, if 
the member belongs to payload PI, then { x pjJ “ { x pl} • 


OF POOR QUA Li's Y 


80 



Therefore, once the system response is known, one can substitute the 
displacement vector {xp^ into Equation (216) and obtain the member 
loads. This direct approach is called the 'Misplacement M method. In 
many ca^es, this is a perfectly valid approach especially if all the 
modes in Equations (32) and (39) are kept. However, if modes are 
truncated according to a cut-off frequency, this procedure often leads to 
inaccurate results. This can be corrected by using the so-called 
"acceleration" method whereby £ xp \ is replaced in terms of applied 
forces and accelerations using the system equations of motion without 
damping. 


Hruda and Jones introduced a load transformation consistent with 
modal synthesis techniques. In terms of the present notation, 

Equation (216) can be replaced by: 


’■?i\ 


{ V 


= [LTl] 


J -I 






with 


+ [LT2] {x”j} 


am. - t^ j tVpW J il -v (i p$ J ! v 


(217) 


(218) 


[LT2]. = [ttWjlTpjl 


(219) 


Note that [ LT2 ] = 0 when the interface is determinate. 

Normally, the displacement vector ^ w in Equation (217) must be 
written in terms of forces and accelerations using the second partition 
in Equation (193). This not only increases the ## cost of the procedure, 
but also introduces the booster accelerations ^q^} and forces •£ Fg| 
into the Problem. We shall now show that this is not necessary under 
certain conditions. 

First, let us obtain the expression for which ordinarily 

should be used in Equation (217). From the second partition in the 
discrete Equation (193) we obtain for -[x^ : 

Ki- tv 1 «v T w- - < 2 *» 

- tpilsi - > 


0F POOR 
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where 


[B] - [tJm b I b ] 
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I?1 


"piVpi 1 ° 

I 

0 

I 


(222) 


This is the exact expression for the discrete interface displacement 
vector. The damping force [DjI^ £x^ possibly can be neglected, but 
we may also keep it. Next, the acceleration approach requires the 
introduction of the modal accelerations which then subsequently are 
truncated: 


{*?}- lK n rl < IT « ]I { p ,l - - 


- t«n) {^iV ' V 4^1) 


(223) 


This is the expression which ordinarily must be used in Equation (217), 
i.e., the {xfj: as given by Equation (223) must be used to select out 

the appropriate { xfj \ partition. 

On the other hand, the expression for £x^^ as calculated by the 
response routine can be obtained from the second partition of 
Equation (196): 

f*ii ■ ly " 1 < [T B ,T { r Blr - - 

- » n i tff) 

which turns out to oe identical to Equation (223)* Therefore, we can 
use { x Ijl i° Equation (217) as given by the response routine. No 
conversion is necessary. 

In case interface modes are used in the response routine, i.e. if 




( 225 ) 
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then the calculated vector is given by: 

Ki- Kr l ( rt>“l T tT B l T {F B J. 
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The corresponding -£xj|is given by: 

{A\ - 1$ r 1 ( [4 *j] t it b ] t {f b } 

- - t2J,W 1 l 

- n B j) 


(226) 


(227) 


where we used Equation (225). 

Expression (227) will be identical to expression (223) if and only if 
all interface mode3 are kept. Indeed , considering the properties 


^®] T [M II ]^] = M J , ^] x [K IX ]^] 


B,T, 




we can write 


[K 


II 


l" 1 = $\) 1 1 


(228) 


(229) 


only if we keep all interlace modes. Conversely, if we do not keep all 
interface modes, then Equation (229) is not valid and the right-hand side 
of Equation (22 ^ is different from the right-hand side of 
Equation (226). Substituting Equations (225), (228) and (229) into 
Equation (227) yields: 


l*\l- tv T Kf- 




B-,-1 


(230) 
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or, 
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{«'}- [Kjj ]- 1 < H b J T { F b } - - 

- m{?J- t«nir*;h V(^) 

which is identical to Equation (223) 

Next, let us write Equation (218) as follows: 

[LTl]^ * [k^tlpj] - [EpjlllpjMpjIpjll^ ] 
fE Pj 1 f I Pj M Pj T Pj 1 ^ 


(231) 


(232) 


It is now relatively easy to show that [E p ^ ] [ I^M p jI p j) ] = 

p; -1 

even when not all the payload modes are retained in 

# N 

P 

[^J. Therefore* Equation (234) can be written as: 

[LTl], = [kfljilpjJ ( - t^l t 


[EpjlUpjMpjTpj] ) 


(233) 


Using Equations (21 ^ , (219) and (233) has several advantages over 
the approach outlined in Reference£53l * First* in case the interface is 
indeterminate (i.e. if LT2 ■ 0) it is not necessary to use Equation (223) 
for | x® £ , if and only if we keep all the interface modes [^] . 

This not only makes the evaluation of Equation (217) much simpler, but 
also reduces the amount of information to be stored in the course of the 
response calculations (i.e* only^q^j-, end £x^£ must b * 

stored). Secondly, the rather .expensive evaluation of (of]* * 
-wi j ii pj iiE !>j ni- j M p .i pj ]^i .. proposed in Reference ££*3.] 

is now replaced by the more efficient computation of l f3 )j * 

$>o£j] 1 . Indeed, applying unit loads successveiy to each of 

the non-interface dofs., often becomes a rather expensive item, 


considering the potentially very large payload mcdtls. The other term in 
Equation (232), (T]j = j t Ipj ) (E p j ] * 

xfl^jMpjTpj] can be easily evaluated by first forming the product 
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[ipjMpjTpj], the columns of which can be looked upon as inertial 
loads applied at the interface of the payload* The corresponding 
deflections are equal to [£ ] . ■ [E„ . ] [ iT .M^ .T„ . ] • Finally, one 
must evaluate [£ ]j = [T^] - n Equation (219) so that we obtain: 

I LTl ) j « t <T] j) (234) 

[LT2] - £]j (235) 

Observe that the payload organization can easily save 
[ <5* ] j , and [£ ]; so t^at any member in ''he payload can now be 
investigated without recalculating these uu_ ; ces. Observe that 
Reference[53}requires the payload analyst to make the choice of members 
before the evaluation of Equations (234) and (235). If for some reason 
an additional member has to be investigated, a reevaluation of [ LTl ] j 
and [LT2]j is necessary. 

Finally, it should also be noted that the loads calculation does not 
involve the "modal modes'* which reduces the computational cost even more. 

5 . Conclusions and Summary 

A numerical integration scheme has been presented. It is a 
"full-scale" approach in the sense that it does not introduce new 
assumptions or approximations compared to the convent ional "exact" solution 
techniques. Improvements over the conventional techniques are introduced 
in both the response and loads calculations. 

The response analysis uses an adaptation of the Newmark-Chan-Beta 
numerical integration technique. This integration scheme is directly 
applied to the coupled system equations (i.e. booster/payload system) 
thereby avoiding the expensive solution of a system eigenvalue problem. 

The Newmark- ‘Chan-Bet a technique has the convenient feature that the step 
size can based on the "cut-off frequency" associated with the forcing 
function regardless of the highest system f requency. This particular 
feature is necr ^sary in the present method because the highest system 
frequency is not known a priori. Although there are techniques to 
determine the highest frequency, it is very likely that this highest 
frequency will be much larger than the cut-off frequency which would lead 
to a much smaller time ste^. It should also be noted that the present 
approach allows for the solution of much larger systems. 

Next, we derived a load transformation consistent with the above 
modal synthesis method. Several cost saving features vere introduced. 
First, we showed that in the case of an indeterminate interface it is not 
necessary to write the interface displacements in terms of accelerations 
and forces, provided one keeps all the interface modes. As indicated, 
one could actually keep the discrete interface displacements instead of 
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introducing the interface modes. Secondly, it was shown that several 
simplifications can be affected in the calculation of [LT1] and [LT2] 
leading to a more efficient and convenient loads calculation. Finally, 
it should also be pointed out that wc do not involve system modes which 
reduces the cost and simplifies the analysis. 

It is estimated that the present approach will reduce the computer 
cost of a payload integration effort by a considerable amount. 

Considering the numerous load cases that must be considered ; n the course 
of a design effort, the present approach may prove to be of ;raat value. 
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CHAPTER III. A PAYLOAD INTEGRATION TECHNIQUE - 
SHORT-CUT VERSION 


1. Introduction 


In Chapter II, we developed a numerical integration scheme which is 
capable of solving the set of coupled system equations as given by 
Equations (196). The method is a full-scale technique because no special 
assumptions or approximations are used. This means that the results are 
as accurate as those obtained from a standard full-scale approach. 
However, under the right circumstances, tr.e cost of this approach is 
significantly lower. 

In this chapter, we wish to adapt this full-scale direct integration 
technique to a short-cut version. The approach is based on estimating 
the size of the feedback from the payload response into the booster. 
First, we shall derive the so-called coupled base motion equations for 
the system. These equations still represent a set of accurate full-scale 
equations of motion. One possible approach is to completely neglect the 
feedback of the payload into the booster. This leads to a technique 
known as the direct base drive technique as discussed in Chapter I, 
section A3f and B3. In many cases this technique leaas to acceptable 
results. Instead of completely neglecting the feedback, we thill subject 
the magnitude of this feedback to a criterion. The result will be a 
method which, depending on the nature of the structure at hand, will 
vasilate between a full-up ^oapled base motion technique and a direct 
base drive method. 

In section 2^ we shall derive rhe pertinent equations of motion. 

2 . The Equations of Motion 

First , note that the equations of motion for the coupled base moti 
technique are derived in Chapter A3*f. Following that derivation, let us 
write down the coupled set of equations of met ion as given by 
Equations (196): 
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This set of equations represent the equations of motion for the 
booster without payioad(s). Note that one could also load the booster 
with a standard payload which would represent a kind of average payload. 

The solution of Equation (238) is the same for all pa/load 
configurations an' does not change as long as the forcing function and 
tiie booster model do not change. A modified Newmark-Chan-Beta numerical 
integration scheme can be used to obtain the solution of (238): 


{iloi+i * t’ioi + (1 - 7)h «t 01 + »"hlVf 0i+1 
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Pr emul t iply ing (240) by [ -BD^ I I] yields 
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The equations to use then a. Equations (239, 241, 243, 244). 
quantities to save are » £ x lof ^*10^ 
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Returning to Equation (236) we can write: 
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g . Q 

The residual quantities {q^Rr *nd £*irJ are clearly due to the 
presence of the payload(s), i.e. they are the feedback of the payload(s) 
back into the booster. Let us write the following vector equation. 
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Substituting Equation (246) into Equation (236) and taking into 
account Equations (237-238) leads to the coupled base motion equations, 
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where 
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(248) 


This set of Equations (247) could be solved resulting in a full-scale 
accurate solution. However, physically it is possible that the feedback 
vector £ sma H (note that we need not have and 

( X Ir£ small), i.e. 0 for all times t. Then from the third 

partition of Equation (247) we obtain a decoupled equation for <[q^£ , 


♦t«ra 
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which can be easily solved, because \ x iof i* Renown. This approach is 
called the direct base drive technique and has been used successfully. 

The problem however is that the magnitude of J xjrj is not known 
in advance and may not always be small. In that case. Equation (247) 
rhould be solved retaining the coupling terms. Again, a Newmark-Chcn- 
Beta technique can be used as follows, 


i+i " + ( 1 _y)h V < ij 1 + ^ h i^} 1+ i 

+ + (0.5-£)h 2 \q} i + h2 ft{h’} 1+1 


(250) 
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( 255 ) 


At this point, it is possible to introduce a criterion which checks 
the magnitude of say • If this magnitude is smaller than a 

certain preset £ then the quantities A J{x® i+1 } and A j£x^ i+1 ^ 

in Equations ( 2 53 ) are not calculated. A possible criterion could be 
of the following form : 



( 256 ) 


where £ is a preset percentage (e.g. 0.01). 
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This criterion could partially avoid the premultiplications by A£ 
and a£ in Equation (253)* There are two more cost generating 
premultiplications by A 2 and A 3 in Equation (252)* These could 
possibly also be avoided when the feedback acceleration is small. This 
would mean a direct base drive at that particular time step. Both these 
approaches were implemented and will be discussed in the next chapter. 

The main problem with this kind of approaches is to find an answer to 
the question: What constitutes a small feedback? This question is still 

not answered even with an equation like Equation (256). Even though 
encouraging results were obtained, it is recognised that additional 
research and development is necessary. 
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CHAPTER IV. THE SOFTWARE PACKAGE - IMPLEMENTATION 


1. Introduction 


This chapter discusses in general terms the software package 
associated with a complete booster /payload response and loads analysis. 

An attempt will be made to clearly link the theory of Chapter I with the 
specific program and subroutine descriptions. This will give us the 
opportunity to touch upon some of the constraints and difficulties 
invariably associated with the development of a practical payload 
integration software package. Some factors to consider are: computer 

core usage; convergence; available data; the separation of booster, 
payload and integration organizations; work schedules; engineering time; 
ease of program usage; computer cost and related efficiency of 
algorithms; reuse of existing information; required accuracy versus cost; 
handling of potentially large models ; etc. 

Section 2 of this chapter presents a general description of the 
organization and components of the software package. In particular, we 
explain the purpose and contents of the components and how they relate to 
each other. 

Section 3 presents a simple sample problem and shows how it is 
analyzed and evaluated. Also, the case of the STS-ST-OMS Kit is 
discussed in addition to some other sample problems. Finally, we shall 
also discuss some of the results related to the short-cut version 
developed In Chapter III. 

2. Organization - General Description 

This section outlines the organization of the software package. 

Figure 4 represents a flow diagram of a complete booster/payload 
integration problem. Each of the flow diagram blocks has a program 
associated with it. Therefore, there are six programs: PROGRAM BOOSTER, 

PROGRAM PAYLOAD, PROGRAM INTFACE, PROGRAM FORCE, PROGRAM RESPONS AND 
PROGRAM LOADS. Each of these programs draws on a pool of subroutines 
called FORMA (Fortran Matrix Analysis). FORMA is a library of 
subroutines coded in FORTRAN IV for the efficient solution of structural 
dynamics problems. These subroutines are in the form of building blocks 
that can be put together to solve a large variety of structural dynamics 
problems. The FORMA library was developed by the Dynamics and Loads 
Section at Martin Marietta Aerospace and is being updated and expanded 
whenever the need occurs. 

It should be pointed out that other libraries can be used and that 
the proposed integration method does in no way inherently depend on the 
FORMA library. However, in this report, the software is built around the 
FORMA subroutine library (In particular, the Partition-Logic version) and 
therefore, the * r is assumed to have a working knowledge of that 
library. 
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Figure 4: PAYLOAD INTEGRATION PROBLEM 
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There are several reasons motivating the PROGRAM approach. Because 
all FORMA routines are written in terms of variable dimensions, it is 
possible to write each PROGRAM for the specific dimensions of the problem 
at hand, thereby optimizing computer core usage. Also, the user often 
has at his/her disposal d-’ta already generated by other means. The 
PROGRAM approach allows l..± user to omit the recalculation of that 
particular data. For example, a set of so-called expanded modes could be 
available. The user then can directly read those expanded modes into the 
PROGRAM and omit the use of a subroutine which calculates those expanded 
modes. Furthermore, the PROGRAM approach allows for the separation of 
data generated by the booster, payload and integration organizations. 
Indeed, very often these three organizations are physically at different 
locations and data peculiar to one organization often is not readily 
available to the other organizations. PROGRAM BOOSTER for example, only 
deals with data pertaining to the booster and, therefore, works 
independent of the payload. Also, sometimes the data generated by 
PROGRAM BOOSTER can be used in analyses of different payloads and 
therefore has not to be recalculated. 

Finally, the PROGRAM approach also allows for better check-out and 
control of the data generated at several points in the process of a load 
cycle. Indeed, the user can put in his/her own checks if desired. 

PROGRAM RESPONS represents the hub around which the other five 
programs are centered. The purpose of PROGRAM RESPONS is to generate the 
coupled booster/payload system response. The most important subroutine 
called by PROGRAM RESPONS is SUBROUTINE ZRESP, which implements the 
integration sc’ ae as outlined in Chapter II. The INPUT to PROGRAM 
RESPONS consists of all the quantities necessary to run SUBROUTINE 
ZRESP. The OUTPUT of PROGRAM RESPONS is the system response i.e. 
displacements, velocities and accelerations. These quantities can be 
written on paper and tape. In particular, the payload accelerations and 
the booster/payload interface displacements at each time step are written 
on tape so that they can be used in PROGRAM LOADS for the calculation of 
member loads. 

Much of the INPUT to PROGRAM RESPONS is not directly given and 
therefore must be created in advance. PROGRAM BOOSTER, PROGRAM PAYLOAD 
and PROGRAM FORCE were composed to serve this purpose. PROGRAM BOOSTER 
generates all booster data necessary to run PROGRAM RESPONSE. The 
booster organization car use this PROGRAM independently of any other 
organization. Enough subroutines were developed so that all booster data 
can be generated stirtlng with the free mass and stiffness matrices 
[Mol and [Kg] and the interface restrained modes and frequencies 
[^f] and [ £•§]. It is reasonable to expect that these INPUT 
quantities are available. If not, the user is expected to provide this 
information before running PROGRAM BOOSTER. It woul^ not be wise to 
"can” the construction of [Mg], [Kg], [<^^1 and [<3g) because of 
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the multitude of ways these quantities can be generated. Furthermore, 
PROGRAM BOOSTER contains a number of "flags” which allow for user 
flexibility of INPUT. Indeed, often times certain quantities are already 
available and need not to be regenerated. The same is true for PROGRAM 
PAYLOAD which is very similar to PROGRAM BOOSTER except that it generates 
payload quantities necessary to run PROGRAM RESPONS. In addition, it 
also generates parts of load transformations if desired. Again, much 
flexibility is possible depending on the case at hand. PROGRAM INTFACE 
collects some of the data generated by PROGRAM BOOSTER and PROGRAM 
PAYLOAD and produces quantities that involve both booster and payload 
data. Again, these quantities are needed in PROGRAM RESPONS and PROGRAM 
LOADS. PROGRAM INTFACE reflects the coupling between booster and payload 
through the Interface. For example, ic calculates the interface modes 
[$l]. PROGRAM FORCE essentially converts the force data into the 
right format to be used in the integration program PROGRAM RESPONS. 
PROGRAM FORCE also contains a number of "flags” which allows for more 
flexibility. Finally, as mentioned above, PROGRAM LOADS generates member 
loads and draws on PROGRAM PAYLOAD for load transformation INPUT and on 
PROGRAM RESPONS for payload response INPUT. 

Each of the six PROGRAMS are independent components of the software 
package. PROGRAM BOOSTER can be used by an Independent booster 
organization. Similarly, PROGRAM PAYLOAD can be used by an independent 
payload organization. PROt AM INTFACE, PROGRAM FORCE and PROGRAM RESPONS 
can be used by an independent integration organization while PROGRAM 
LOADS can be used by any organization that is responsible for loads 
calculations. Because each of the PROGRAMS is compatible with the other 
PROGRAMS, it is also possible for one organization to use the entire 
package in sequence. 

This section was intended to give the reader a general idea of how 
the software package is structured. It is not intended to be a detailed 
user guide. The Final Report of this contract will be accompanied by a 
detailed user guide as well as the actual listings of all the PROGRAMS 
and associated SUBROUTINES. 

3. Numerical Examples 


In this sectirn we shall discuss several simple sample problems which 
were used to check out the internal correctness of the software package. 
Furthermore, we shall briefly present the results of two realistic 
analyses namely, the analysis of the STS-ST-OMS Kit system and the 
analysis of a defense booster/payload system. 

The first example Is depicted in Figure 5. The booster B consists of 
18 pipe segments. The mass of each segment is equally divided between 
the end points of the segment. If we only keep translational dofs, then 
the free booster h;»s 57 dofs and the "cantilevered” booster has NB - 54 
dofs. Similarly, the payload P consists of 7 pipe segments and NP - 21. 
Because there are 3 rigid body modes, we have a determinate interface and 
IF - 3. The parameters for a booster pipe segment are: 

E-6.89xl0 6 kN/m 2 , J-2.77xx0 2 kg/m 3 ,A body -1.93xl0- 2 m 2 , 

A^Hg-4. 17x10“ ^m 2 ,L“0. 762m J 0 «8. 325xl0“Sn^ 
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Similarly, for a payload pipe segment: 

E*6. 89xl0^kN/m 2 , -2. 77xl0 2 kg/m 3 ,A= 4.05xlO*V, 

L“0. 762m, J 0 «*l. 249x1 0~6 m 4 

Using the above data, a finite element model was derived for both the 
cantilevered booster and the cantilevered payload. Solving the 
eigenvalue problem yields booster frequencies ranging from 1 Hz to 106 Hz 
and payload frequencies from 1 Hz to 104 Hz, In this particular example 
we used zero Initial conditions and applied loads to stations 16 to 17 in 
the x, y and z-directions (444822xcos(i50t)N,i*l,6) . 

The accuracy of the adapted Newmark-Chan-Beta routine was checked by 
comparing the response results from r OGRAM RESPONS with those obtained 
from a fourth order Runge-Kutta (Gili modification) routine using the 
same step size. The results compare very well. Table I shows some of 
the results for a step size h ™ 0.001 seconds. It should be pointed out 
that a phase shift in the response was observed. This is to be expected 
and is inherent to the numerical technique used. However, this phase 
shift does not affect the value of the maximum or minimum loads In an 
element. It may slightly affect the time at which this minimum occurs, 
but this is of little consequence. Indeed, it is fair to say that in 
practice every numerical integration technique produces a phase shift 
when a time step is used consistent with a cut-off frequency. 

Next, we compared computer cost of the present method with that of 
the conventional approach. The conventional "full-scale" approach first 
calculates the so-called "modal modes and frequencies" (i.e. the system 
modes from Equation(196) ) , Then, a numerical scheme (e.g. Runge-Kutta) 
is used to determine the response from the uncoupled system equations. 

For the present example, the cost of the direct integration routine to 
determine the response is less by a factor of 5 compared to the 
conventional approach. 

We also compared the cost of the load calculations. The improved 
technique decreases the cost by a factor of 6, It Is hard to tell how 
this factor will change when the booster/payload system represents a more 
re ‘stic configuration. Also, it may be hard to compare cost factors 
for large systems from a logistic point of view. Indeed, payload 
organizations currently deliver load transformations which are not 
consistent with the present improved approach. Therefore, it is often 
impossible to generate the appropriate quantities required for use in the 
present approach, for lack of certain information. It should be noted 
that this is only a logistics problem and changes could easily be 
accommodated. 

Table II lists some of the load results obtained for this example. 
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Th* next temple problem contltte of e booster plate model to be 
integrated with two truss-like payload*. The finite element code FINBL 
(a finite element program developed by the Martin Marietta Aerospace 
Dynamics Sectlou) was used to develop the three models for the booster 
and the two payloads. 


Th* booster as ihown in Figure 6 is a hexagonal cylinder 
conaistlng of 24 quadrilateral plat* sections. Th* material properties 
for these plates are: 

9 - 0.025 lb/ln 2 , K - 10.6 x 10* lb/in 2 
V ■ 0.334 , t (thickness) ■ 0.1 In 


Th* 24 elements were joined at 30 nodes as described In Table 1. 
Bach of these 30 nodes were assigned three trsnslatlonal degrees of 
freedom. Th* geometry and dof. numbering scheme for this booster model 
are shown In Table 2. Twelve nodes (numbers <3, 9, 11, 12, 14, 15, 17, 
18, 20, 21, 23 and 24) and all of their correapondlng dofs. were 
designated to make up the interface for possible coupling with payloads. 
Thus, for '"he booster, we have 36 Interface dofs. (IF ■ 36) and 54 non- 
interface dofs. (Nd ' S4). The booster model was forced in the 
x-dlr*ction at nodes 26, 27, 29 *nd 30 (dofs. 78, 81, 87 and 90) as 
follows: 


l.e. 


r 78 

f 82 

f 87 

f 90 

NF 


150000 * SIN(27T (0.5) (t-0.001) ) 
160000 * SIN[ 27T (0.45) (t-0.0u2)J 
125000 * SIN[27T(0.6) (t-0.C2)) 
17000* SIN(2TT (0.55) (t-0.01)l 
4. 


Figure 7 show* the first payload PI, which is made up of 18 oar 
elements. The bars are joined at 8 nodes to fora this truss. The 
geometry of the structure is given in Table 4 and the material properties 
in Table 5. All nodes are assigned translational degrees of freedom 
only. The four comer nodes make up th* interface. We have: ND * 12, 

IF ■ 12, NP • 12 and th* Interface dofs. are 1, 2, 3, 4, 5, 6, 10, 11, 

12, 13, 14 and 15. 


The model of payload P2 is shown in Figure 8 . It consists of 28 bars 

4 olned at 12 nodes forming a truss. Th* geometry and th* degree of 
freedom table for the model are given in Table 7 and the material 
properties and connections in Table 8. Again, only translational dofs. 
are considered. The parameters are: ND - 36, IF • 12, NP « 24 and the 
interface dofs. are 7, 8, 9, 10, 11, 12, 19, 20, 21, 22, 23 and 24. 

A total description of how to use tne software package for this 
sample problem is given in the User Guide MCR-82-&02 
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Figure 6 , Booster model for Sample Problem 2. 
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INPUT OAT A FOB CON8INEO MEMBRANE -BENDING QUADS I LATERAL PLATE ELEMENTS 


MASS 

RO 

T(MASS) 


Ml STIF - K1 

. 250E *01 £ * .106E+06 

. 100E+00 NU * . 334E+00 

T( MEMBRANE ) » - 100E+00 

T ( BENDING ) - . 100E+00 


ELEMENT 

JOINT 1 

JOINT 2 

JOINT 3 

JOINT 4 

NUMBER 





i 

1 

2 

8 

7 

2 

2 

3 

9 

6 

3 

3 

4 

10 

9 

4 

4 

5 

11 

10 

5 

5 

6 

12 

11 

6 

6 

1 

7 

12 

? 

7 

8 

14 

13 

8 

8 

9 

15 

14 

9 

9 

10 

16 

15 

10 

10 

1 1 

17 

16 

1 1 

1 1 

12 

18 

17 

12 

12 

7 

13 

18 

13 

13 

14 

20 

19 

14 

14 

15 

21 

20 

15 

15 

16 

22 

21 

16 

16 

17 

23 

22 

17 

17 

18 

24 

23 

18 

1b 

13 

19 

24 

19 

19 

20 

26 

25 

20 

20 

21 

27 

26 

21 

21 

22 

28 

27 

22 

22 

23 

29 

28 

23 

23 

24 

30 

29 

24 

24 

19 

25 

30 


Table 1 , Description of quadrilateral plate 
elements for the booster model used 
in sample problem 2 
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JOINT DAT* USED IN SUBROUTINE FENKA 



DEGREES OF 
TRANSLATION 

FREEDOM 

rotation 


GLOBAL CARTESIAN COORDINATES 

JOINT 

U 

V 

V 

P 

0 

R 

X 

y 

Z 

1 

1 

2 

3 

0 

0 

0 

5.0000 

0.0000 

50.0000 

2 

4 

5 

6 

0 

0 

0 

2.5000 

4.3300 

50.0000 

3 

7 

8 

9 

0 

0 

0 

-2.5000 

4.3300 

50.0000 

4 

10 

1 1 

12 

0 

0 

0 

-5.0000 

0.0000 

50.0000 

5 

13 

14 

15 

0 

0 

0 

-2.5000 

-4.3300 

50.0000 

6 

16 

17 

18 

0 

0 

0 

2.5000 

-4.3300 

50.0000 

7 

19 

20 

21 

0 

0 

0 

5.0000 

0.0000 

25.0000 

8 

22 

23 

24 

0 

0 

0 

2.5000 

4.3300 

25.0000 

9 

25 

26 

27 

o 

0 

0 

-2.5000 

4.3300 

25.0000 

10 

28 

29 

30 

o 

0 

0 

-5.0000 

0.0000 

25.0000 

1 1 

3 1 

32 

33 

0 

0 

0 

-2.5000 

-4.3300 

25.0000 

12 

34 

35 

36 

0 

0 

0 

2.5000 

-4.3300 

25.0000 

13 

37 

38 

39 

0 

0 

0 

5.0000 

0.0000 

0.0000 

14 

40 

41 

42 

0 

0 

0 

2.5000 

4.3300 

0.0000 

15 

43 

44 

45 

0 

0 

0 

-2.5000 

4.3300 

0.0000 

16 

46 

47 

48 

0 

0 

0 

-5.0000 

0.0000 

o.ooco 

17 

49 

50 

51 

0 

0 

0 

-2.5000 

-4.3300 

0.0000 

18 

52 

53 

54 

0 

0 

0 

2.5000 

-4.3300 

0.0000 

19 

55 

56 

57 

0 

0 

0 

5.0000 

0.0000 

-25.0000 

20 

58 

59 

60 

0 

0 

0 

2.5000 

4.3300 

-25.0000 

21 

61 

62 

63 

0 

0 

0 

-2.5000 

4.3300 

-25.0000 

22 

64 

65 

66 

0 

0 

0 

-5.0000 

0.0000 

-25.0000 

23 

67 

68 

69 

0 

0 

0 

-2.5000 

-4.3300 

-25.0000 

24 

70 

71 

72 

0 

0 

0 

2.5000 

-4.3300 

-25.0000 

25 

73 

74 

75 

0 

0 

0 

5.0000 

0.0000 

-50.0000 

26 

76 

77 

78 

0 

0 

0 

2.5000 

4.3300 

-50.0000 

27 

79 

80 

81 

0 

0 

0 

-2.5000 

4.3300 

-50.0000 

28 

82 

83 

84 

0 

0 

0 

-5.0000 

0.0000 

-50.0000 

29 

85 

86 

87 

0 

0 

0 

-2.5000 

-4.3300 

-50.0000 

30 

88 

89 

90 

0 

0 

0 

2.5000 

-4.3300 

-50.0000 


Table 

2. 

Geometry description and degree of freedom table 
for the booster model used in sample problem CL 
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ALL BAR ELEMENTS 
p = 0.0025 lb/in 3 
E = 10.6 x 10 3 lb/in 2 
G = 3.84 x 10 6 lb/in 2 
A = 0.01 in 2 
J Q = 1.67 x 10" 5 in 4 
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Figure 7 . Payload 1 



JOINT data used in subroutine femka 


oR 't'S 

Of * 


i V 




DEGREES 

C c 

FREEDOM 



global 

CARTESIAN coordinates 


TRANSLATION 

RC* 

AT ION 





JOINT 

u 

V 

w 

P 

n 

R 

X 

f 

z 

1 

1 

2 

3 

C 

0 

0 

2.5000 

4.3300 

0-0000 

2 

4 

5 

€ 

0 

0 

0 

- 2 5000 

4.3300 

0 0000 

3 

7 

8 

O 

0 

0 

0 

•2.5000 

- 1 oooo 

1 .0000 

4 

10 

1 1 

* 

' <■ 

0 

0 

0 

-2.5000 

-4. 3300 

c.0000 

5 

13 

14 

IS 

n 

0 

0 

2.5000 

-4 3300 

e.ocoo 

€ 

16 

17 

18 

Q 

0 

0 

2 . 5000 

• 1 . 0000 

1 .oooo 

7 

19 

20 

21 

o 

0 

0 

-2.5000 

- 1 .0000 

- 1 ocoo 

8 

22 

23 

24 

0 

0 

0 

2. 5000 

-1.0000 

- 1 -OOOO 


Table 4* Geometry and degree of freedom table 
for payload l 


INPUT DATA POR BAR ELEMENTS 

KODEK * KDDEB = 

MASS * Ml STIF = LOAD TRANS = PAV 3L T STRESS TRANS = 

RO * .250E-02 E = . 106E+05 ALPHA * o. 

G = 394E+07 


element 

JOINT 

JOINT 

REF 

area 

POLAR 

TORSION 

z bending 

> BENDING 

SHEAR 

NUMBER 

1 

2 

POINT 


inertia 

const 

INERTIA 

inertia 

factor 

1 

6 

1 

2 

1 OOE - 0 4 

167E-04 

. 4 4 1 E -04 

.B33E-C5 

833E-05 

.333 

2 

1 

2 

3 

1 0'"' E - 2 1 

. 167E-04 

. *4 "04 

S33E-05 

.833E-05 

.833 

3 

2 

3 

1 

. i OOE -01 

. 167E-04 

. <4!E-04 

833E-05 

333E-05 

.833 

4 

3 

6 

1 

. 1 0OE *01 

’67E-04 

. 4 4 1 E - 04 

833E-05 

. S33E-OS 

.833 

5 

3 

1 

2 

100 E -0 1 

'676*04 

. 141E-04 

.933E-05 

.833E-05 

.833 

6 

8 

2 

i 

. 1 OOE -0* 

167E-04 

4 4 IE -04 

833E-05 

.833E-05 

.833 

7 

3 

4 

€ 

. 100E-E i 

. 167E-04 

. 141E *04 

.833E-05 

833E-05 

.833 

8 

4 

5 

i 

lOOE-O 4 

* 67 E -Q4 

. 14 IE *04 

.833E-05 

. 8 3 3 E - 05 

.833 

9 

5 

6 

3 

1 OCE -0 1 

157E-04 

. * 4 1 E '04 

.833E-05 

.833E-C5 

833 

io 

4 

6 

1 

10 CE * 0 1 

167E-04 

. # 41E~04 

833E-05 

.833E-05 

.833 

1 1 

7 

5 

4 

. 1 OOE - 0 1 

. 167 E *04 

. 14 IE *04 

.833E-05 

.833E-05 

.833 

12 

6 

8 

5 

. 1OCE-01 

. 167E-04 

. 14 IE *04 

833E-05 

833E-05 

.833 

13 

3 

7 

4 

lOOE-Oi 

. 167E-04 

. 141E-04 

833E-05 

833E ^5 

.833 

14 

8 

7 

2 

i OOE '01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

15 

5 

8 

€ 

. 100 E - 0 « 

. 1G~E-04 

1 <dM E *04 

.833E-05 

.833E-05 

.833 

16 

4 

7 

3 

100E-01 

. 167E-04 

. 141E-04 

833E-05 

833E-05 

.833 

17 

1 

8 

6 

100 E -01 

. 167E-04 

. 4 4 1 E -04 

.933E-05 

.833E-05 

. 833 

18 

2 

7 

3 

100E -0" 

. 167E-04 

14 1 t -04 

.833E-05 

.833E-05 

833 


Table 5* Payload 1 material properties 
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JOINT DATA USED IN SUBROUTINE FEMKA 


orsyv" is 

Jr f\ > .JAUTY 


JOINT 


OEGREES OF FREEOOM 
TRANSLATION ROTATION 

U V W P 0 R 


GLOBAL CARTESIAN COORDINATES 
X Y Z 


1 

1 

2 

3 


4 

5 

6 

3 

7 

8 

9 

4 

10 

1 1 

12 

5 

13 

14 

15 

6 

16 

17 

18 

7 

19 

20 

21 

8 

22 

23 

24 

9 

25 

26 

27 

10 

28 

29 

30 

1 1 

31 

32 

33 

12 

34 

35 

36 


OOO 
0 0 0 

0 0 0 

0 0 0 

OOO 
0 0 0 

OOO 
OOO 
OOO 
OOO 
OOO 
OOO 


2.5000 

2.5000 

2.5000 

-2.5000 

-2.5000 

-2.5000 

-2.5000 

2.5000 

-2.5000 

2.5000 

-2.5000 

2.5000 


0.0000 

2.0000 
4.3300 
4.3300 

2.0000 
0 0000 
-4.3300 
-4.3300 
0.0000 
0.0000 

2.0000 

2.0000 


1.0000 
1 .0000 
0-0000 
0.0000 

1.0000 

1.0000 
0-0000 
o.oooo 
- 1.0000 
- 1.0000 

- 1 oooo 

- 1 .0000 


Table 7. Geometry and degree of freedom table 
for payload 2 


INPUT DATA FOR BAR ELEMENTS 




KODEK 

* 

KODEB - 

MASS - Ml 

STIF * K 1 

LOAD TRANS * 

PAY4LT 

STRESS TRANS * 

RO * 

.250E-02 E « 

G « 

. 106E+05 
. 384E+07 

ALPHA * 

0 . 


ELEMENT 

JOINT 

JOINT 

REF 

AREA 

POLAR 

TORSION 

Z BENDING 

Y BENDING 

SHEAR 

NUMBER 

1 

2 

POINT 


INERT I.* 1 

CONST 

INERTIA 

INERTIA 

FACTOR 

1 

1 

2 

4 

. 100E-01 

. 167E-04 

. 14 1 E *04 

.833E-05 

.833E-05 

.833 

2 

2 

3 

4 

. 100E-01 

. 167E-04 

, 14 1 E -04 

.833E-05 

.833E-05 

.833 

3 

3 

4 

2 

. 100E-01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

4 

4 

5 

2 

. 100E-01 

. 167E-04 

. 14 IE-04 

.833E-05 

.833E-05 

.833 

5 

5 

6 

1 

. 100E-01 

. 167E-04 

. 14 IE-04 

.833E-05 

.833E-05 

.833 

6 

6 

7 

1 

. 100E-01 

. 167E-04 

. 1 4 1 E -04 

.833E-05 

.833E-05 

.833 

7 

7 

8 

1 

. 100E-01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

8 

8 

1 

6 

. 1005-01 

. 167E-04 

. 14 IE -04 

.833E-05 

.833E-05 

.833 

9 

5 

2 

3 

. 100E-01 

. 167E-04 

. 14 IE -04 

833E-05 

.833E-05 

.833 

10 

1 1 

12 

10 

. 100E-C1 

. 167E-04 

. 1 4 1 E “04 

.833E-05 

.833E-05 

.833 

1 1 

5 

3 

2 

. 100E-01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

12 

6 

1 

2 

. 100C-01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

13 

6 

2 

1 

100E-01 

. 167E-04 

. 1 % IE-04 

833E-05 

.833E-05 

.833 

14 

9 

10 

12 

. 100E-01 

. 167E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

15 

7 

1 

8 

. 100E-01 

. 167E-04 

14 IE -04 

.833E-05 

.833E-05 

.833 

16 

4 

12 

3 

. 100F-01 

. 167E-04 

. 141E-04 

.833E-05 

833E-05 

.833 

17 

1 1 

10 

9 

. 100E-01 

. 167E-04 

. 14 IE -04 

833E-05 

833E-05 

.833 

18 

9 

8 

10 

. 100E-01 

. 167E-04 

. 14 IE -04 

G33E-05 

833E-05 

.833 

19 

1 

10 

12 

. 1O0E-01 

. I67E-04 

. 141E-04 

.833E-05 

.833E-05 

.833 

20 

2 

12 

10 

. 100E-01 

. 167E-04 

. 14 IE -04 

.833E-05 

.833E-05 

.633 

21 

6 

i 

1 1 

. 100E-01 

. 1 6 7 E - 04 

. 14 IE -04 

.833E-05 

.833E-05 

.833 

22 

5 

1 1 

9 

. 1C0E-01 

. 167c -04 

. 141E-04 

.833E-05 

. 0336 -05 

.833 

23 

8 

10 

1 

. 1O0E-01 

. 157E-04 

. 141E-04 

.833E-05 

833E-05 

.833 

24 

10 

12 

1 

. 100E-01 

. 167E-04 

. 14 IE-04 

.033E-O5 

.833E-05 

.833 

25 

12 

3 

2 

. 100E-0I 

. 167E-04 

. 14 IE-04 

.833E-05 

.833E-05 

.833 

26 

7 

9 

6 

. 100E-01 

. 167E-04 

. 14 1 E -04 

833E-05 

.833E-05 

.833 

27 

J 

1 1 

6 

. 100E-01 

. 167E-04 

. 14 IE-04 

.833E-05 

.833E-05 

833 

28 

1 1 

4 

5 

. 100E-01 

. 167E-04 

. 14 IE -04 

833E-05 

.833E-05 

.833 


Table 8. Payload 2 material properties 
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The integration was carried out over 0.9 seconds* i.e. STARTT * 0, 

ENDT * 0.9 and DELTAT - 0.01. 

The results are encouraging. The cost per time step is equal to 
0.00018 C.U. (C.U. ■ cost unit) in subroutine ZRESP whereas the cost per 
time step for the usual TRSP3 is 0.00017 C.U. It should be recognized, 
of course, that the usual TRSP3 route requires the solution of a 126 x 
126 eigenvalue problem which costs 0.031 C.U. The accuracy of the 
response is the same for both approaches. The loads program ZLOADS 
requires on the average 4.x 10"® C.U. in order to track a maximum and 
minimum load at a particular station. This compares with 9.X10"^ C.U. 
for the conventional technique. Again, the maxima and minima compare 
well with the values generated by the conventional technique. 

The theory as presented in Chapter III was also applied to the case 
of the MX missile. Indeed, the method lends itself very well to this 
case. The missile vas basically modeled in two parts, an aft and front 
■>nd. The number of interface dofs. is 6 being the ideal number for this 
integration technique. Also, a respectable number of load cases must be 
evaluated. In particular, several payload configurations have to be 
investigated. The conventional technique requires a new system 
eigenvalue problem for every new configuration of the Reentry Vehicles. 

This undertaking was conceived as a parallel effort to the regular 
techniques. Again, the results were encouraging. It should be noted 
that damping was included, resulting in slightly different results for 
the loads (differenced 0.2%). This can be attributed to the difference 
in handling the modal damping. The model consisted of a 261 degree of free- 
dom payload model and a 92 degree of freedom booster model. A cut-off 
frequency of 50 Hz was chosen, resulting in a coupled booster/payload 
model containing 27 cantilevered booster modes and 55 cantilevered payload 
modes. The entire analysis was performed on the VAX/VMS-1 1/780 computer 
sytem. 


Another realistic sample problem is given by the Space Transportation 
System/Space Telescope/OMS Kit structure. Figure 9 represents the 
Space Telescope. 

The payload for the Space Telescope (S.T.) mission consists of two 
cargo elements: (1) the Space Telescope and (2) the OMS Kit. These two 

structures are coupled to the Space Transportation System (S.T.S.). 
Together, these three structures form the lift off system. The S.T. 
itself is a modal synthesis of the System Support Module (S.S.M.) and the 
Orbital Telescope Assembly (O.T.A.) using interface restrained modes, A 
free S.T. mass and stiffness matrix of size 214 x 214 was formed and 175 
cantilevered mode shapes and frequencies calculated. There are 6 
interface degrees of freedom with the S.T.S. This means that the 
interface is determinate. 

The OMS Kit is modeled with 36 degrees of freedom and has 7 degrees 
of freedom in common with the S.T.S. 
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The S.T.S. lift off dynamic model for the hot S.R.B. propellant 
condition is again provided in terms of interface restrained 
coordinates. The model has 759 degrees of freedom and 300 interface 
restrained modes and frequencies were calculated. The S.T.S. interface 
has 30 degrees of freedom. This means there are 30-7-6 “ 20 superfluous 
interface coordinates on the booster side. 

A modal synthesis of the lift off system (S.T.S. + S.T. + OHS Kit) 
was accomplished through a Craig-Bampton formulation and is used for 
comparison purposes. Finally, also provided was a lift off forcing 
function (L07201) which turns out to be a critical case for the Space 
Telescope. 

At this point all necessary data are available to run the six program 
software package. 

The number of booster restrained modes kept in the analys* 

NB “ 300. Similarly for payloads 1 and 2 we have NP1 175,. 36. 

The interface dimensions are IF « 30, IF1 * 6, TF2 =» 7, and Z . 0 . 

The number of non-zero applied forces on the booster is equal 1. 

The possible cost savings occur in programs RESPONS and LOADS. For 
the current S.T.S.-S.T.-0MS Kit system we considered a start time equal 
to zero and an end time equal to 10 seconds. The integration step is 
equal to 0.005 seconds. Furthermore, the velocities and the elastic 
accelerations at time zero are chosen to be zero. 

A CDC/CYBER 75C was used to execute all programs. In this particular 
case, program RESPONS required on the average 0.0015 C.U. (“realistic 
cost unit) per time step h. It should be noted that a phase shift was 
observed in the response compared to the conventional approach where 
system modes ("model modes") are used. This is typical for a 
Newmark-Chan-Beta numerical scheme. This phenomenon is no cause for 
concern if the only goal is to arrive at maximum and minimum member 
loads. Although it is hard to correctly compare different techniques, 
the authors feel that the cost to obtain the response compares favorably 
even with the conventional system modes approach where a set of decoupled 
system equations are solved. The average cost per time step here becomes 
0.0017 C.U.’s. Note that when system modes are used, one has to perform 
an additional [C*>] T type multiplication on the right hand-side of the 
decoupled equations. In this case this results in a mere expensive force 
term calculation than is the case in the present scheme. 

The maximum and minimum member loads obtained from program LOADS are 
all within 1% of the corresponding loads obtained from the conventional 
approach. The gain in program LOADS is primarily due to the fact that no 
system modes and frequencies are involved in the computations. In 
addition, the term [LT2] j can be evaluated directly without 

having to write £xjjJ- in terms of accelerations and applied forces as 
should be done when system modes are used. The average cost in order to 
track the maximum and minimum load at a particular station Is 0.00012 
C.U. which compares with 0.00026 C.U. 1 8 for the convention?.! technique. 


114 



Finally, we wish to discuss some of the results obtained for the 
short-cut versions as discussed in Chapter IV. We tried out two schemes 
based on Equation (256) in Chapter IV. The actual control statement used 
can be written as follows: 



where £ is a small quantity chosen by the user. The left-hand side of 
the inequality represents a normalized magnitude squared of the interface 
acceleration feedback. Tne first scheme checks this quantity against the 
chosen £. If the inequality Is fulfilled, then the calculations of 
the A 2 and A 3 t^rms In Equation (253) are by-passed. The second 
scheme follows the philosophy of a direct base drive approach. This 
means that if the inequal Is fulfilled at time t, we assume that the 
feedback is small at that time, and a direct base drive solution is 
generated for that time. If the Inequality is not fulfilled, a coupled 
solution is obtained for that time step. That means that this short-cut 
versic vasilates between a full-scale and a direct base drive method. 

The first scheme turned out to be expensive. The reason for this is 
that a coupled base motion technqiue requires a rather significant amount 

of bookkeeping . 

The second scheme proved successful in the case where the feedback is 
Indeed small. Again, because of the added bookkeeping, it is necessary 
that the £-inequality is indeed fulfilled part of the time, because 
otherwise the method becomes more expensive. It is also 1 allzed that 
more research and checking is necessary before a final judgement can be 
given. Compared to the ZPESP routine, the second short-cut version 
vasilates anywhere between 40% slower (coupling at all times) and 30% 
faster (no coupling at all). 

It should be noted that ZRESP (fully coupled) is already within the 
realm of reasonable cost and additional cost reduction, although 
possible, may complicate matters too much. 
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CHAPTER V: CONCLUSIONS AND SUGGESTIONS 


I. Introduction 


In Chapter I-A we reviewed and assessed a set of "full-scale” 
methodologies. This allowed us to introduce the necessary background 
material in terms of a unified nomenclature. All these methods have 
their merits. However, the Residi:*;! Mass and Stiffness Method appears to 
yield the most accurate results. It is the approach which best describes 
the booster structure in terms of a minimum number of modes, given a 
certain cut-off frequency for the externally applied force £ Fg | • The 
fact that no payload Information is required to obtain the booster model 
is a very convenient feature in connection with the present study. 
Therefore, the same booster model can be used as long as the booster does 
not change. 

In C* .pter I-B several short-cut methods have been discussed and 
evaluated. Although each of these methods has its own merits. It is 
believed that none of the is acceptable in their present stage of 
development to function as a standard short-cut method for general use. 

In fact, It would be very hard to develop an "ultimate" short-cut 
method. The main reason for this may be the multitude of different 
situations such a method would have to accommodate . 

In Chapter II we attempted the development of a full-scale numerical 
integration scheme to obtain the response of a booster/payload(s) 
system. It is felt that under the right circumstances this method can be 
of certain value. 

A short cut version of the above numerical scheme was developed in 
Chapter III. The method Is designed to vasllate between a full-up 

upled base motion approach and a direct base drive method. It Is hoped 
that this technique will have the characteristics of both methods, i.e., 
the accuracy of the coupled base motion technique and the speed of a 
direct base drive. 

In Chapter IV we discuss the implementation and programming of the 
above ujethods. A software package was developed and checked out on 
several sample problems. 

In the remainder of this last chapter we shall discuss a few other 
suggestions for possible short-cut approaches. 
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2 . Bass Not too Technique* 


ORIGINAL l - 
OF POOR QUALHY 


The Coupled Base Notion Technique as explained In Chapter I-A, leads 
to a fundamental aet of equations, which we repeat here for clarity of 
presentat ton, 
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As mentioned In Chapter l-A, the payload designer Is prlstarlly 
Interested In the response of the payload l,e,. 
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where £Sf^(and {A^^must be computed from Equations ,1 57-259) and 
t rom Kquatton (90). The base drive method tocuses o .quatlon (258) 
whtch vtelds provided on the right hand side of Eouatlon 

(258) is known. The tdea Is tt> produce an expression for }M**l wi thout 
actually ao Ivins the couplod aet of Kquatton* (2*7-259). * 

Ft rat, consider the coefficient a o f {A*} and In Equation 

(25*)), These coefficients represent the ratio of the payload mass and 
the total vehicle mass, In many STS applications tM* ratio will be 
rather small « 10X). Therefore, a first possibility la to Ignore these 
terms tn Equation (259), Secondly, In many appllcat tons we can assume a 
statically determinate interface, 1 * * • f Tm KiJm J " m C0J , sc 

that Equation (259) becomes PO 
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Ordinarily* the eoefflctent matrix of Equation (262) is ngt. 

small and cannot he Ignored A first possibility is to assume that 

is small and can he ignored, Thia means that the feedback of the payload 
is not Important. This can be a realistic assumption because the payload 
is usually small compared to the booster. In this case we can completely 
ignore in Equation (258) and write 


Equation (263) is now effectively decoupled fr, m Equations (257) and 
(259). Physically, ignoring the feedback of the payload means that 
payload and booster are not modally coupled. Equation (263) is the 
direct base drive equation. 


A second possibility is to scale the vector j in Equation (262). 
Indeed, let us assume a full-scale solution is available for aome payload 
Pj. Now, some relatively small changes are made in the payload Pj to 
generate payload P. The assumption now, is that $**} is uot much 
different from 

^ X 
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or 


Equation (262) for payload P, can be written as 

- - [’kv'jivt.] “ [*,-.] 
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Taking into account Equation (264) It fc! *..*»• from Equations (266) 
and (262) that 
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which yield* a acaled value for JX^. j to be us*d In Equation (258)* 
Again, one should Investigate when such an approach !* valid. 
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3. Another Possible Approach 


ORIGINAL PAGE IS 

OF POOR QUALITY 

In the course of our investigation and evaluation ot several 

short-cut methods it was noted that many methods involve assumptions and 
approximations leading to either doubtful or cumbersome results* In 
addition, it is often very difficult to assess the effects of those 
assumptions on the response and the loads of the boos ter /pay load system. 


The basic problem Is to somehow deal with the couplir*g effects 
between booster B and payload P without solving an eigenvalue problem 
pertaining to the coupled booster/payload system* This is a difficult 
problem indeed* Each of the short-cut methods discussed in Chapter I-B 
addresses this problem in a different way. However, the proposed 
solutions invariably lead to cumbersome mathematics and program coding* 
This observation led us to the development of a more direct approach 
which we think shows great promise* This new approach Is easy to 
understand and easy to Implement* It Is based on the work of C. W. White 
and B. D. Maytum [74], 


Let us recall Equation (25) l 
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which represents the set of equations of motion of the coupled 
booster/payload systwm. It is now assumed that a cut-off frequency is 
defined based on a Fourier series expansion of Furthermore, we 

also assume that e.g. the Residual Mass and Stiffness Method was used to 
construct the following set of modally coupled equations* 
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and the cut-off frequency was used to letermine the size of [ ] and 

J. In other words, the size of Equation (269) is already much less 
than the size of Equation (268). Due to e.g. the Residual Mass and 
Stiffness Method, the reduced Equation (269) still represents an 
acceptable model for the coupled booster /pa/load system. 


The first step of the present approach is to solve the eigenvalue 
problem associated with Equation (269), namely 
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yielding a set of modes [ ] and a set of frequencies 

satisfying 
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and 
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where ^uj* are the new normal coordinates. Substituting transformation 
(272) into Equation (269) and premultiplying by [4*]^ and using 
Equations (273-274) we obtain the uncoupled s^t of equations, 



(275) 


The modal matrix [ty ] and the frequency matrix represent the 

modal information of the coupled booster/payload system. The idea now is 
to change the payload and calculate the changes in [(f'] and f/l?J . In 
other words, we use the full-scale solution of Equation (269) as a 
"start-solution'*. This approach is taken in most short-cut methods and 
as such does not detract from the present approach. For example, this 
full-scale solution could be determined at the beginning of a design 
effort and would stay the same for all subsequent design cycles of a 
particular payload. 
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OF POOR QUALITY 

Let us now consider the new payload Pj., with mass matrix [KpjJ 
and stiffness matrix [MpiJ. This payload P^ could be totally new or 
just a modification of the nominal payload P, as long as we have the same 
degrees of freedom for both payloads P and Pj. For this new payload 
Pj, we replace Equation (268). 

, 1/ .. i r v ?B1 \ 

^!!?i^i!sijipi!Ipi!pi_ UJ + K H^ r ?i K pi T Pi j 0 ( _* 2 i_ ^ (276) 
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is the new transformation matrix and 
_*B1_ 

"Ipi" (278) 

*N 

is the new system displacement vector. 

Let us no«v write Equation (276) as follow^: 
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( 282 ) 


( 283 ) 


(234) 


Note that in case the interface is statically determinate 
[k BB ] “ [0] and in addition the geometry of the payload is not 
changed then [TpjJ ■ [Tp ] and 
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Also, note that if no changes are made in the mass the right hand 
sides of Equations (282) and (285-288) become [0] and similarly, if no 
stiffness changes are made we have from Equation (277) that 
[Tp] ■ [Tpj] and consequently Equations (285-288) are valid while 
Equation (289) becomes [kpp] * [0] although the interface can still be 
statically indeterminate. 


Next, let us define the following transformation 
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After substituting Equation (290) into Equation (279) and 
premultiplying by [A]* we obtain the set of equations that now replaces 
the set Equation (269) 
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The next step is to define the transformation 
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The transformation (292) is now substituted i ...o Equation (291) after 
which we premultiply by [vf'] T and invoke properties (273-274), yielding 
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This Equation (293) replaces Equation (275). For convenience, let us 
denote 
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Matrices [M] and [kl represent the perturbations in the mass and 
stiffness matrices fIJ and of system (275). At this point, 

several observations can be made* First it should be noted that it is 
very possible that certain changes in the payload will only affect a 
limited number of modes and frequencies. This means that several columns 
in [V] and corresponding elements in will not change after the 

changes in the payload are made. This reduces the size of Equation 
(293). Secondly, in solving the eigenvalue problem associated with 
Equation (293) it is possible to use a Rayleigh-Ritz approach with fIJ as 
the estimated start modes. The smaller the changes in the payload the 
better estimate [Ij will be and the less iterations will be necessary to 
produce the new modes and frequencies of the perturbed booster B/Payload 
PI system. An even better starting set of modes could be the solution to 
the perturbed eigenvalue problem with all off-diagonal terms equal to 
zero (this is equivalent to the first term in a Taylor series expansion 
of the perturbed system modes and frequency). Thirdly, we wish to 
investigate the possibility of truncating modes in [*¥ ] according to the 
initially defined cut-off frequency. If this was possible Equation (293) 
could be reduced in size by approximately 50% compared to the already 
reduced system Equation (269). This reduction would be in addition to 
the one due to unaffected modes as mentioned above. However, this 
question must still be carefully investigated. Finally, It is also 
possible that the modes are grouped in subsets which show very little or 
no coupling between each other. This means that the eigenvalue problem 
associated w^th Equation (293) can be replaced by two or more smaller 
eigenvalue problems, which of course reduces the computation time. 

There are additional advantages to this method: simplicity of use; 

accuracy of results (e.g., tnis method could even be used as a full-scale 
method); possibility of using engineering judgement and experience; the 
possibility to identify changes requited to meet certain frequency 
requirements; the possibility to change brarichf requencies to decouple the 
load problem leading to smaller eigenvalue problems, the potential for 
significant computational time savings. 
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